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Summary 


This is the user's manual for UCAP (Unified Counter-rotation Aero- 
acoustics Program), the counter-rotation derivative of the UAAP 
(Unified Aero-Acoustic Program). The purpose of this program is to 
predict steady and unsteady air loading on the blades and the noise 
produced by a counter-rotation Prop-Fan. The aerodynamic method is 
based on linear potential theory with corrections for non-linearity 
associated with axial flux induction, vortex lift on th e blades, and 
rotor-to-rotor interference. The theory for acoustics and the theory 
for individual blade loading and wakes are derived in Unified 
Aeroacoustics Analysis for High Speed Turboprop Aerodynamics and 
Noise, Volume 1 (NASA CR-4329) . This user's manual also includes a 
brief explanation of the theory used for the modelling of counter- 
rotation. 


Manor sections of this volume include a general code overview, 
descriptions of input and output, installation instructions, and a 
listing of error codes. 


I Introduction 

This document describes the use and general structure of the 
Unified Counter-rotation propeller Aeroacoustic Program (UCAP) , 
developed under NASA Contract NAS3-24222, Task Order 10. UCAP is 
the single and counter-rotation successor to the single rotation 
UAAP,’ which is a computer program to predict the aerodynamic and 
acoustic performance of Prop-Fans using helicoidal lifting 
surface theory and a frequency-domain acoustic theory developed 
at Hamilton Standard in the 1980's (reference 1). 

The major improvements of this version of UCAP over UAAP are: 

The inclusion of counter-rotation. The steady interaction 
is modelled by perturbing the solution for one rotor by the 
mean velocity field caused by the other. These interaction 
velocities are assumed to follow streamlines. 

Calculation of induced drag from blade Ac p data instead of 
in the Trefftz plane. This was changed to include loading 
from the leading edge vortex and the tip edge vortex. 

Inclusion of the vortex loading effects into the calcula- 
tions of the induced velocity and the wake. The time aver— 
age, non-linear axial momentum is satisfied by iteration. 

Streamline contraction is added. This is modelled by forc- 
ing continuity along annuli defined by flow-field data. 

Steady aerodynamic boundary conditions on the blades are 
calculated on conical surfaces defined by the local stream- 
line angle. 

The calculation of the required turning angle (for the flow 
tangency boundary correlation on the blades) has been re- 
vised. Instead of treating the axisymmetric disturbances 
(for example, those from the other rotor or from an external 
flow field) as perturbation angles, these are incorporated 
into the mean velocity triangle (see Figure 1) . 

The axisymmetric interaction component of the rotor-to-rotor 
interaction velocity is calculated from dCJdx, the local 
thrust loading, and dc p /dx, the local power loading. This 
method requires less CPU time than the near wake formula- 
tions and will not significantly change interaction veloc- 
ities at normal rotor spacings. 


Unified Aero-Acoustic Program. 
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Iterative correction for interaction velocity on both rotors 
of a counter-rotation Prop-Fan is included. 
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II Theory 

UCAP is a derivative of the earlier UAAP work; the details of 
the theory behind UAAP can be found in reference 1. The two major 
improvements over the interim version of UCAP, produced under 
Tasks I and II of this contract, are the incorporation of vortex 
loading into the induced flow calculation and the wake calcula- 
tions, for both single and counter-rotation operation, and the 
modeling of contraction in rotor-rotor interference, which will 
be described below. UCAP works on the basis of perturbing each 
rotor with the velocity field produced by the other rotor. The 
steady loading upon a blade is divided into several parts: the 

potential part (subscript pot) is directly from lifting surface 
theory; it is calculated from 

{Lp OT > = [K) 1 x{W} 

where {L^} is a vector representing the potential loading, [K]' 1 
is the influence coefficient matrix, and {W} is the turning 
angles required to maintain flow tangency along the blade sur- 
face. The secondary loading (that portion of loading due to 
induced flow) is calculated based on momentum theory using the 
thrust and torque at each control point and the mass flux through 
the rotor. The methods for calculating the secondary loading and 
vortex loading (the portion of the loading from the leading edge, 
side edge, and tip edge vortices) are described in reference 1, 
but the application is different. In UAAP, the calculation of 
the induced flow and the wake depended only upon the potential 
loading, {L^} . The steady potential loading was the result of 
non-linear iterations, where the C p and C, resulting from each 
step were used to compute a new induced flow. This was repeated 
until the change in C p was small. The final potential load was 
then used to drive the Near Wake Calculation. The vortex loading 
was attached to the final performance values but never included 
in either the loading to drive the Near Wake Calculation or in 
the induced flow calculation. In the present program, the 
induced flow is calculated before and after the calculation of 
the vortex loads, so the influence of non-linear lift is included 
in the inflow field. The Near Wake routine calculations are 
driven by {L}, which includes the vortex load and is obtained by: 

{L} = {L pot } . ( OCt/ax) /(dcjdx) ) 

where dC^/dx includes the effect of vortex loading and dcjdx 
includes potential loads only. 

The flow chart in Figure 2 shows a schematic of the method used 
in counter-rotation operation for the steady aerodynamic solu- 
tion. Single rotation operation bypasses steps M through V, but 
is otherwise similar. 
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Streamline contraction is a result of the conservation of mass. 

In UCAP this is done by dividing the flow into streamtubes which 
are defined by the flow field entered by means of the VELGRADS 
command (see Section IV, "Input Description"). If no external 
flow field is defined, initial streamlines will be based on the 
hub and tip radii of the front and rear rotors. The locations 
where these streamlines intersect the rotors are calculated after 
every steady iteration. The new locations of the streamlines at 
the rear rotor are computed by calculating the flow rate in each 
streamtube and adjusting the annular areas in order to enforce 
incompressible continuity. For the streamline locations on the 
rear rotor, the flow rate in each streamtube where it crosses the 
forward rotor is calculated from the induced velocity, the 
interference velocity caused by the rear rotor, the velocity due 
to the external flow field, and the streamtube area. Since the 
velocity where that stream-tube crosses the rear rotor is known, 
continuity will determine the stream tube area at that location. 
This is repeated until all streamlines on the rear rotor are 
relocated. The new locations of the streamlines where they 
intersect the front rotor are calculated in a similar manner: the 
streamtube areas on the rear rotor are held fixed while the 
streamlines on the front rotor are adjusted to satisfy continu- 
ity. 

The streamtubes defined by the last steady iteration are used for 
the calculation of the unsteady rotor— to— rotor interferences. A 
separate document (reference 4) contains a more detailed explana- 
tion of the theory used in UCAP. 
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Ill Program Organization 

UCAP is a large program, with about 30000 lines of source code 
excluding comments. It is divided into several modules, which 
communicate via common blocks, argument lists, or files. These 
are: 

Control Module — Controls the program flow, monitors conver- 
gence of the steady aerodynamic iterations, and performs 
some housekeeping functions, such as assuring that the 
common blocks are properly loaded with default data. 

Input Module — Reads the data entered by the user. Entered 
data are placed into named common blocks for access by the 
rest of the program. The input module will also set the 
default values for the appropriate input parameters which 
have them. - ~ 

Blade Geometry Module — Converts the geometry data entered 
in any of the three available coordinate systems to the 
coordinate system used internally. 

Steady Aerodynamics Module — Performs the calculations re- 
quired for the solution of the steady aerodynamics, and 
passes the Ac p 's calculated to the Steady Loads Module, 
which performs the integrations needed to get C , C,, and 
applies the vortex loading to get C* and C*. The influence 
coefficient matrices and thickness vectors are calculated 
here. The forward rotor influence coefficient matrix is 
written to file FWD000. The rear rotor is written to file 
AFT000 . For single rotation, the file will be SRP000. 

Steady Loads Module — Integrates Ac p 's to get C p , C t , dC p /dx, 
and dCJdx, and applies the vortex loading corrections, P as 
in UAAP , to get Cf, C*, dC*/dx, and dc*/dx. This module 
also writes the steady portion of the input data file for 
the Noise Module. 

Unsteady Aerodynamic Module — Calculates the complex Ac 's 
resulting from a given disturbance field and writes the P 
unsteady portion of the input data file for the Noise Mod- 
ule. The unsteady influence coefficient matrices calculated 
in this module are written to external files with names of 
the form AFT001, AFT002, FWD001 , FWD002 , etc., where the 
number is the interaction harmonic order. For single rota- 
tion all matrices will be written to a file named SRP000. 

Steady Interaction Wake Module — Uses dc*/dx and dc*/dx to 
calculate the velocity field caused by one rotor upon the 
other. This approach discards influence of the chordwise 
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loading distribution on the velocity field, which is consid- 
ered negligible at reasonable rotor-to-rotor spacings. 

= '; Near Wake Module — Uses complex loading information supplied 

~ by the Steady Aerodynamics Module, with a correction for the 

vortex loads, to calculate the unsteady perturbation terms 
caused by one rotor upon the other. While the perturbations 
are unsteady, the wake is steady when viewed by an observer 
moving with the rotor generating the wake. When an inflow 
field incorporating unsteady components is input, the UCAP 
wake is based on the steady component only. The unsteady 
portion does not impact the wake. For single rotor opera- 
tion, only, this module will calculate the velocities caused 
r- by the rotor at specified upstream or downstream points 

m (WAKEEXEC command) . 

Noise Module — Calculates the radiated noise from data con- 
^ tained in an input data file at locations requested by the 

user. 

— Utility Routines — Includes sub-programs for interpolation, 

splining, data management, and printout. These routines are 
called from many places in the program. 

w Figure 3A shows the subroutine calling tree of the program. 

Figure 3B shows the subroutine cross reference listing. Figure 
3C is a subroutine reference and purpose listing. Figure 3D 
shows the labeled common area reference listing. 

UCAP uses a large number of named common blocks, but no unnamed 
(blank) common blocks, for data storage and for data exchange 
*■' between modules. These are listed in Table I. Further informa- 

tion concerning the actual system requirements is in Section VII, 
"Program Installation". 
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IV Input Description 

The primary input mode is expected to be a predefined input data 
set. Upon program execution the primary input file, file number 
5, is immediately read and echoed to both the primary output 
file, file number 6, (to obtain a complete copy of the input on 
the output file) and to a scratch file, file number 11. All 
subsequent primary input is obtained by reading file number 11. 

No other input files are required by the code except as the user 
may define for restart capabilities. Several scratch files are 
generated by the program; they are defined in the Section VIII, 
"File Requirements". 

There are many input values which act as tolerances on itera- 
tions, or convergence criteria. These have been defaulted within 
the code to "recommended" values. Although they explicitly 
appear in the input description, they need not be input. Input 
data which need not be input appears at the end of each section, 
and the input location number appears within parentheses. 

The next sections present a description of all the input neces- 
sary to run the code. 
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Commands 


The code is command driven. Commands are used to identify the 
input sections, the input values, and to execute various func- 
tions of the code. Sub-commands are used to execute options 
within functions. The general input format of the command is: 

(starting in column 1) 

COMMANDX ( SUBCOMMD ) where "COMMANDX" is an 8 character word, in- 
cluding trailing blanks, and "SUBCOMMD" is 
also an 8 character word, including trailing 
blanks . 


Load 


Some commands require that numerical input must be read in next. 
For this purpose a location-specified, free-field input routine 
is used which reads between columns 1 and 72. Locations for input 
are indicated by "L" followed by a number indicating the desired 
start of a location. 

Some of the location specified input controls integration mesh 
sizes, Fourier series convergence, and other program tolerances. 
These locations have been identified by parenthesis () around the 
location number and can usually be ignored by the user. However, 
if these values are input, the parenthesis () should be omitted 
from the location field. 

Input for this "LOAD" format is illustrated below: 

(starting in column 1) 

C THIS IS AN OPTIONAL COMMENT RECORD (or Records) 

L 10 1.1 2 3.5 

C THIS IS ANOTHER OPTIONAL COMMENT RECORD (or Records) 

12.7 13.1 

L 91 6.5 

END 

This will cause the locations shown below to have the following 
values: 


3CATION 

VALUE 

10 

1.1 

11 

2.0 

12 

3.5 

13 

12.7 

14 

13.1 

91 

6.5 


Notes: "Scientific notation" is not allowed, e.g. 2.1 E+03 is 

not allowed. Implied repetitions are not allowed, e.g. 3*2.1 is 
not allowed. The "END" record is required to terminate each 
entry into the LOAD routine. 
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Command Summary 


A list of accepted commands is shown below. Since this is a 
command driven code, command order is important, and therefore 
the commands are listed below in the required logical sequence 
for execution. 


HEADER 

RUNPARMS 

CRPPARMS 

AIRPARMS 

BLADEGEO 

LSTPARMS 

NOIZPARM 

VELGRADS 

VORTPARM 

WAKEPARM 

INTERPRM 

AEROEXEC 

WAKEEXEC 

NOIZEXEC 

ENDCASE 

ENDJOB 


Input page header cards 

Input flight parameters 

Input counter-rotation parameters 

Input options for 2-D drag look-up. 

Input propeller /blade geometry (2-D,RXY or 
XYZ coordinates) . 

Input options to panel aero code. 

Input options to the noise portion of the 
code. 

Input axial velocity field on a defined grid. 

Input options to vortex flow calculation. 

Input options to wake calculation. 

Input options for rotor-to-rotor interaction 
calculations 

Execute the aero code. 

Execute the wake calculation. 

Execute the noise portion of the code. 

End of input and calculations for the 
current case. 

End of job, terminate the program. 
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Specific Input Description 

Input requirements /options are provided below for each command, 
in the order listed above. 


HEADER : 

This command enters page heading records. Records following the 
HEADER command are sequentially read until an •'END” command is 
found. Up to 10 page heading records may be entered. 


RUNPARMS : 

This command enters run parameters using ’’LOAD" format. The 
location, default, and input descriptions are: 


LOC. 

Default 

Variable 

1 

0 

RUNBUG 

2 

59.0 

QDEGF 

3 

1.0 

QRHOR 

4 

0 

QADVF 


5 

0 

QMX 

16 

0 

QADVR 


Description 


Debug Option, 1 turns on. 

Ambient temperature, °F. 

Ambient density/ Sea Level 
Std. density. 

Advance ratio, J * V 0 /nD, 
for a single rotor or for 
the forward rotor of a 
counter-rotation Prop- 
Fan. 

Free-stream Mach number. 

Advance ratio, J = V 0 /nD, 
for the aft rotor of a 
counter-rotation Prop- 
Fan. 


CRPPARMS : 

This command is used to control counter-rotation operation of the 
program and to enter certain parameters specific to 
counter-rotation, such as the rotor-to-rotor spacing. If this 
command is omitted, UCAP will run in single rotation mode; this 
is to permit existing UAAP input data to be used for single 
rotation without change. 
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LOC 

1 


2 


3 


4 


6 


Default Variable 

0 . 0 CRPBUG 


0 . 0 SWITCH 


-1.0 SPACE 


0 . 0 COUNT 


1000.0 CRPTOL 


Descriptions 
Printout control. 

0: Minimal printout 

1: Printout specific 

data after each front/r- 
ear rotor iteration 


0.0: Run single ro- 

tation. This was chosen 
to ensure compatibility 
with UAAP data decks. 

1.0: Run counter- 

rotation 

Distance between rotor 
pitch change axes, nor- 
malized by forward rotor 
diameter. For counter- 
rotation cases, this must 
be entered by the user; 
if nothing is entered the 
program will stop. 

Number of front-rear ro- 
tor iterations. Maximum 
is 25. 

Tolerance. When the 
change in power coeffi- 
cient between iterations 
is less than this value, 
the steady performance is 
considered converged. 
Values less than l.OxlO' 5 
will be reset to 10000. 
This value was selected 
for convenience in pro- 
gram debugging. The user 
must enter a value less 
than 1. for iteration to 
occur. See Appendix for 
recommended values. 
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0 . 0 FWDHRM 


8 


0 . 0 AFTHRM 


Highest forward rotor 
wake harmonic to use for 
the excitation of the 
rear rotor. 

Highest rear rotor induc- 
tion field harmonic to 
use for the excitation of 
the forward rotor. 


LSTPARMS : 

This command enters input options to the panel aero portion of 
the code using "LOAD” format. This command has two sub-commands: 
AFT and FWD. Use LSTPARMS (FWD ) to enter parameters for the 

forward rotor in counter-rotation cases. Use LSTPARMS (AFT ) 

to enter parameters for the rear rotor in a counter-rotation 
case. Use LSTPARMS, with no sub- commands , for single rotation 
cases. The location, default value and input description are: 


( * See 

Figure 4 for 

pictorial 

description) 

LOC. 

Default 

Variable 

Description 

1 

0 

ZSTBUG 

Debug Option, 1 turns on 

* 2 

10 

QNCP 

Number of chordwise pan- 
els, maximum of 10. 

* 4 

8 

QNSM 

Number of spanwise modes 
(control point radii) . 
Maximum of 10. 

20 

1 

QPART1 

0 generate [K]' 1 and W x . 

1 generate [K]' 1 and rea< 

W x . 


2 read [K]' 1 and generate 

W x . 

3 read in [K]' 1 and W T 
The code requires an in- 
verse kernel matrix [K]‘* 
and thickness vector, W x , 
to obtain the aerodynamic 
loading on the blade. 

The [K] _l matrix and W x 
vector require a signifi- 
cant amount of CPU time, 
and an option to utilize 
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previously generated 
[K]“‘ matrices and W T vec- 
tors has been incorporat- 
ed into the code. 

Note: For counter-rota- 

tion operation, the pro- 
gram requires that the 
[K] -1 matrix be available 
for every QQ and QK (see 
below). If any [K] -1 
matrix is required for a 
given rotor, all [K] _l 
matrices must be calcu- 
lated for that rotor. 

Order of unsteady load- 
ing, harmonic, = 0 for 
steady loading. 

For single rotation oper- 
ation, this parameter 
must be set to a value 
greater than zero for 
unsteady aerodynamics. 

For counter-rotation, 
enter 0, as this parame- 
ter is calculated for the 
unsteady rotor-to-rotor 
interactions . 

Number of circumferential 
modes for unsteady load- 
ing, = 0 for steady load- 
ing. 

For single rotation oper- 
ation, this parameter 
must be set to a value 
greater than zero for 
unsteady aerodynamics. 

For counter-rotation, 
enter 0, as this parame- 
ter is calculated for the 
unsteady rotor-to-rotor 
interactions . 


28 


2 


QNBOPT 


29 
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QITNON 


* 101 .2, .35, QZAR 

.45, .55, 

.65, .75, 

.85, .95 


* 141 .5,. 5, QCONTP 

. 5 , . 5 , 

. 5 , . 5 
.5, .5 


181 


. 4 , . 4 QCHW 

.4, .4 
.4, .4, 

.4, .4, 

• 4 , . 4 , 


= 1 for supersonic lead- 
ing edge element, use 
when supersonic flow ex- 
pected at leading edge. 

= 2 for subsonic leading 
edge element. 

Max number of non-linear 
iterations, for single 
rotation operation. 

=0 for linear cal- 
culation. 

This parameter is ignored 
for counter-rotation op- 
eration. 

Spanwise locations of 
control point 
radii. 

The number of these must 
agree with QNSM, fraction 
of R^. 

Chordwise location of 
control points within 
each panel normalized to 
the panel width. 

There must be QNSM of 
these. All control 
points are at the same 
location at a given radi- 
us. 


Width for chordwise 
averaging of 
downwash, normalized to 
panel width. There must 
be QNSM of these input. 
The downwash averaging 
width is constant at a 
given radius. 


15 



51 0, 0 QWMU Complex downwash vector 

for flowfield for ex- 
ternal to propeller un- 
steady loading calculati- 
on. These can be obtained 
from Section 13, equa- 
tion (13.8) of reference 
1. A routine to generate 
these is in the appendix. 
The values of QWMU must 
be input sequentially 
starting in location 351, 
as real part, imaginary 
part for each control 
point across the 1st 
spanwise station, and 
then proceeding outward 
along the blade span. 
Thus, there must be 
(QNCP*QNSM) pairs of val- 
ues input for QWMU. 


The 

following 

input locations 

may be ignored . 

(3) 

.02 

QDR 

Radial extent of singular 
integral at control point 
radius for wake calcula- 
tion radius. 

(5) 

1024 

QN 

Number of points in FFT 
for terms in kernel inte- 




gration, max. of 2048, 
must be a power of 2 . 

(6) 

.004 

QDELTA 

Axial step size for FFT. 

(7) 

4 

QMODOP 

Spanwise mode shape op- 
tion 

(use default value) . 

(8) 

.0001 

QTOLF 

Tolerance for W (omega) 
integration. 

(9) 

.001 

QTOLT 

Tolerance for other inte- 
grations. 

(10) 

.005 

QTOLS 

Tolerance for summations 

(11) 

20 

QMM1 

Loop limit for harmonic 
summation in wake kernel. 
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(12) 

30 

" " QMM2 

(13) 

10 

QMM3 

(14) 

— 

QITABK 

(15) 


QITABT 


(16) 

(17) 

0 

QINT4 

QPRINT 

(18) 

0 

QPRIN1 

(19) 

0 

QPRIN2 

(23) 

.025 

QKDOWN 

(24) 

.010 

QKSTART 

(25) 

5 

QMM4 

(26) 

.9 

QMBLEN 


Loop limit for harmonic 
summation in bound ker- 
nel. 

Loop limit for harmonic 
summation in thickness 
vector . 

File number allocated for 
[K] _1 matrix storage. 

The default is 8 for sin- 
gle rotation or for the 
FWD sub-command and 18 
for the AFT sub-command. 

File number allocated for 
W T (thickness vector) 
storage. The default is 
9 for single rotation or 
for the FWD sub-command 
and 19 for the AFT sub- 
command . 

Not currently used. 

Additional debug print - 
not recommended. 

Additional debug print - 
not recommended. 

Additional debug print - 
not recommended. 

Radial step size for non- 
-singular bound kernel 

Width of singular region 
for bound kernel integra- 
tion 

Loop limit for summation 
in sound power calcula- 
tion. 

Trailing edge effective 
Mach number for blending 
to supersonic trailing 
edge elements. 
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(27) 

1024 

QNO 

Number of points in FFT 
integration for n = 0 
term in kernel, 
max « 2048, must be a 
power of 2 . 

(221) 

.002 

QINMES 

Step size for radial in- 
tegration in wake kernel . 


BLADEGEO ; 

This command enters the propeller/blade geometry using the "LOAD" 
format. For any propeller and flight condition, the shaft power 
required, and the thrust produced by the propeller are a function 
of the blade angle. In using this code it is recommended that the 
blade angle be adjusted so that the power or thrust calculated by 
this program matches a desired value, (see Sections 3 and 6 of 
reference 2 for further explanation) . 

Nine sub-commands exist to allow input of the geometry in differ- 
ent forms: 

For single rotation these are: 

BLADEGEO ( 2 -DCOORD ) , BLADEGEO (RXYCOORD) , and BLADEGEO (XYZCOORD) . 
For the forward rotor in counter-rotation these are: 

BLADEGEO ( 2 -DCOFWD ) , BLADEGEO (RXYCOFWD) , and BLADEGEO (XYZCOFWD) . 
For the aft rotor in counter-rotation these are: 

BLADEGEO ( 2 -DCOAFT ) , BLADEGEO (RXYCOAFT) , and BLADEGEO (XYZCOAFT) 

The 2 -DCOORD, RXYCOORD, and XYZCOORD sub-commands must not be 
used in counter-rotation operation. The 2-DCOFWD, RXYCOFWD, 
XYZCOFWD, 2 -DCOAFT, RXYCOAFT, and XYZCOAFT sub-commands must not 
be used in single rotation. These are described below. 

BLADEGEO j 2 -DCOORD ) 

BLADEGEO j 2 -DCOFWD ) 

BLADEGEO j 2 -DCOAFT) : 

These sub-commands all use the same form of the blade geometry, 
which is expected to be the most widely used. The 2 -DCOORD sub- 
command is used for single rotation cases, 2-DCOFWD is used for 
the forward rotor of a counter-rotation Prop-Fan, and 2 -DCOAFT is 
used for the aft rotor of a counter-rotation Prop-Fan. Input 
takes the form of the spanwise variation of thickness ratio, 
chord/ diameter ratio, twist, airfoil section designation and 
stacking axis coordinates. This command (sub-command) will 
calculate the blade surface coordinates and interpolate these 
into a form required by the other parts of the program. To 
insure that there are no errors in the blade output description, 
only a limited amount of extrapolation of the input blade coordi- 
nates is allowed. Thus, it is best to provide input stations, X, 
and streamline angles, SLA, which will define the root and tip 
sections of the blade such that the output stations, ZBLDST , can 
be interpolated and not extrapolated. The code requires exactly 
ten (10) inputs defining X. The input described below is illus- 
trated in Figures 5 and 6. 

Note that there are empty locations in the input. These are not 
used by the program. 
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LOC. Default 


Variable 


Description 


31 

32 

33 

41-50 

51-60 

61-70 

71-80 

81-90 


347 


348 


0 

0 

0 

0 


BLDBUG 

BLADE 

D 

SCO 

X 

HOB 

BOD 


CLD 


DTHET 


THTDES 


THTCUT 


Debug option, 0 is off, 1 
is on. 

Number of blades 

Propeller diameter, ft. 

Propeller hub/tip ratio 

Spanwise input stations, 
10 are required, frac- 
tions of R^p. 

Spanwise airfoil maximum 
thickness/chord ratio. 

Spanwise chord/diameter 
variation. 

Spanwise variation of 
design lift coefficients. 

Spanwise twist variation, 
in degrees. 

Twist should be input 
such that twist at the 
75% radius is 0. 

The 75% radius value of 
blade angle at which the 
blade is to be "de- 
signed", degrees. 

The 75% radius value of 
blade angle at which this 
calculation is to be run, 
degrees. 
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381 


711 


721 


731 


900 


390 0 SLA 


720 0 XSWP 


730 0 YSWP 


740 0 ZSWP 


2 1 ZNPCOV 


Note: Propeller blades 

are assumed to be "de- 
signed" at one value of 
blade angle. The input 
values of thickness/ 
chord, chord/ diameter , 
camber, twist and stack- 
ing line are assumed to 
be defined at the value 
of blade angle in Loca- 
tion 347. Location 348 
defines the blade angle 
at which this calculation 
is to be run. 

Spanwise variation of 
streamline angle, 
degrees. Blade airfoil 
sections are assumed to 
be on cones which approx- 
imate the streamlines 
through the propeller. 
This input is the cone 
half-angle, positive as 
shown in Figure 5. 

Spanwise variation of X 
coordinate of mid-chord 
stacking line, fraction 
of R,^, see Figure 6. 

Spanwise variation of Y 
coordinate of mid-chord 
stacking line, fraction 
of R^, see Figure 6. 

Spanwise variation of Z 
coordinate of mid-chord 
stacking line, fraction 
of R^, see Figure 6. 

Number of output stations 
in the chordwise direc- 
tion, max 49. 
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901-949 


0 . , 5 . , PCTCHD 


10, 

15, 

20 

25, 

30, 

35, 

40, 

45, 

50, 

55, 

60, 

65, 

70, 

75, 

80, 

85, 

90, 

95, 

100 
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ZNIS 


951-999 .10,. 30, ZBLDST 

. 40 , . 50 , 

.55, .60, 

.65, .70, 

.75, .80, 

.825, .85, 

.875, .90, 

.925, .95, 

.96, .97, 

.98, .99 

1.0 

The following input locations may 
(346) 2 SWPOPT 


(349) 1 ZKCUT 


(741-751) 1 BAFL 


Chordwise location of 
output stations. Both 
upper and lower airfoil 
surfaces are output at 
the same chordwise 
locations, % chord. 


Number of output stations 
in the spanwise dir- 
ection, max of 49. 

Spanwise location 
stations at which 
the blade surface will be 
defined, radius ratios, 
r/IV 


be ignored . 

Input sweep option. The 
default value is re- 
quired. 

Type of airfoil section 
defined 

0 planar, 

1 conical, (default) . 

2 cylindrical 

Integer characterizing 
the 2-D airfoil at each 
spanwise station. Only 
NACA 16 series airfoils 
can be generated with 
this deck. 
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BLADEGEO { RXYCOORD^ 

BLADEGEO (RXYCOFWD1 
BLADEGEO fRXYCOAFT) : 

An optional method for input of the blade geometry has been 
provided. This form allows for up to 50 input stations. 
Additionally, it requires a table of the displacement of the mean 
camber line from the chord line. The RXYCOORD sub-command is 
used to enter data for single rotation cases. The RXYCOFWD and 
RXYCOAFT sub-commands are used to enter data for the forward and 
aft rotor, respectively, of a counter-rotation Prop-Fan. This 
option requires that the blade sections being described be on 
cylinders whose axis is the centerline of rotation. The input is 
illustrated in Figures 7 through 9. 


LOC. Default 

Variable 

Description 

1 

0 

BLSBUG 

Debug Option, 0 is off, 1 
is on. 

2 

0 

STANO 

Number of spanwise input 
stations, max of 50. 

3 

0 

PCTNO 

Number of chordwise input 
stations, max of 50. 

4 

0 

DIAMET 

Propeller diameter, ft. 

5 

0 

SPINNER 

Propeller hub/tip ratio. 

6 

0 

BLADES 

Number of propeller 
blades. 

51-100 

0 

CTSTA 

Spanwise input stations, 

r/ IV 

101-150 

0 

PCTCD 

Chordwise input stations, 
% chord. 

At each 

spanwise 

station 


151-200 

0 

THKOB 

Maximum blade thickness, 
fraction of chord. 

201-250 

0 

CHDOD 

Blade chord, fraction of 
diameter. 
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251-300 0 


CAMBR 


301-350 0 TWIST 

351-400 0 XMC 

401-450 0 YMC 

451-500 0 ZMC 

501-550 0 XSLA 


Equivalent NACA Series 16 
camber. Obtained by tak- 
ing the non-d imensional 
maximum height of the 
blade mean camber line / 
.05515 , fraction of 
chord . 

Blade section chord an- 
gle, deg. 

X-coordinate of mid-chord 
stacking line, fraction 
of R^. 

Y-coordinate of mid-chord 
stacking line, fraction 

of IV 

Z-coordinate of mid-chord 
stacking line, fraction 

of IV 

Streamline angle, deg. 
This is the flow angle 
relative to the center- 
line of rotation, and is 
used in the sweep angle 
calculation. 


Blade 

mean camber 

line displacement 

table . 

1003 

0 

PCTND 

Number of chordwise 
points in mean camber 
line displacement table. 

1004 

0 

STAND 

Number of spanwise points 
in mean camber line dis- 
placement table. 

1005 

0 

CAMLN 

Ascending array of chord 


fractions, PCTND values. 
Followed immediately by : 


24 



1005 

+ 

PCTND 

0 

CAMLN 

1005 

+ 

0 

CAMLN (I ,K) 
1=1, STAND 

PCTND+STAND 

K=l, PCTND 


Ascending array of 
spanwise radii, fraction 
of R^. Starting after 
last input location, 

STAND values. 

Followed immediately by : 

Array (PCTND* STAND) of 
mean camber line 
displacements 
from chord line at each 
radial station, at 1st 
chordwise location, frac- 
tion of chord. This is 
followed by a similar 
array for the 2nd chord- 
wise location, and con- 
tinues through the PCTND 
chordwise location. 


BLADEGEO (XYZCOORD) 

BLADEGEO j XYZCOFWD1 
BLADEGEO fXYZCOAFTl : 

This form of blade geometry inputs the blade description in XYZ 
coordinates of each blade surface. The XYZCOORD sub-command is 
to be used for single rotation. For counter-rotation operation 
use the XYZCOFWD sub-command to enter data for the forward rotor; 
use the XYZCOAFT sub-command to enter data for the aft rotor. 
Figure 10 illustrates the coordinate system used. Note that this 
XYZ coordinate system is different than the 2-D or the RXY 
system . 


LOC. 

Default 

Variable 

Description 

1 

0 

BLDBUG 

Print option. 

0 is minimal 

1 is additional. 

2 

0 

DIAME 

Propeller Diameter, ft. 

3 

0 

SPINN 

Hub to tip ratio. 

4 

0 

BLADE 

Number of blades 

5 

0 

STINN 

Number of radial stations 
in input. 
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6 0 CDINN 

7 0 STOUTN 

8 0 CDOUTN 

51 0 STOUTV 

101 0 CDOUTV 


Number of chordwise sta- 
tions in input. 

Number of output radial 
stations needed to define 
the blade. 

Number of output chord- 
wise stations needed to 
define the blade. 

Values of radial sta- 
tions, fraction of radi- 
us. 


Values of chordwise sta- 
tions fraction of chord. 


An "END" record is required to terminate the "LOAD” input. The 
"END" record is immediately followed by the xyz coordinates of 
the blade as follows: 

For each of the STINN radial input stations the following input 
records are required: 

a) A label record 

b) For each of the CDINN chordwise input stations the 
X, Y, and Z coordinates of the face (pressure) side 
of the blade are input. With one set of X, Y, Z, 
coordinates in free-field format per record. The 
units are inches. 

c) Another label record 

d) For each of the CDINN chordwise input stations the 

X, Y, and Z coordinates of the camber (suction) side 
of the blade surface are input. With one set of X, 

Y, Z coordinates in free-field format per record. 

The units are inches. 
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VELGRADS ; 

This command will input or initialize the axial velocity ratio at 
the propeller. These ratio may be due to the spinner-hub and or 
nacelle effects on the freestream flow, but the result must 
appear to the code as an axisymmetric flowfield. 


LOC. Default Variable 


Description 


1 
2 

0 : initialize velocity 
field to 1.0, fraction of 
freestream velocity. No 
further input is 
required. 

1 : input the velocity 
ratio. 


0 

0 


V1BUG 

V10PT 


Debug option 

Velocity ratio input op- 




For option V10PT=1. only: The axisymmetric velocity 

ratio is input on a grid 
of radial and axial 
points which encompass 
the blade outer bound- 
aries. The flowfield 
values are interpolated 
at specific "nodal” 
points required by the 
code. 


3 0 VINRD 

4 0 VINAS 

5 0 VIVRF 

6 0 VITIP 


Number of radial stations 
in grid. 

Number of axial stations 
in grid. 

Reference velocity, by 
which local velocities 
are normalized, same 
units as VIRD, VI W, and 
VIAX. 

Reference radius, by 
which local radii are 
normalized, same units as 
VIRD and VIAX. 
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7 


0 


VICLL 


Centerline location. The 
axial location of the 
forward (or only) blade 
pitch axis is assumed to 
be 0. This value is sub- 
tracted from the input 
values of VIAX after in- 
put so that the axial 
location of the forward 
(or only) blade center- 
line will be at 0. Same 
units as VIAX. 


The following inputs are repeated for each of the axial 
stations in the grid, (K=1,VINAS). 
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0 

VIRD ( I , K) 
1=1 , VINRD 

Radial locations of 1st 
set of up-stream grid 
up-stream grid points, 
same units as VITIP. 

25+ 


VI W (I , K) 
1=1, VINRD 

Velocities at VIAX (I) , 
VIRD ( I ) , 1=1, VINRD same 
units as VIVRF. 

25+ 


VIAX ( I , K) 
1=1, VINRD 

Axial locations of 
VIW(I), same units as 
VITIP. 

WAKEPARM: 

This command provides input parameters to the propeller wake 
calculation procedure. This input is used in both the perfor- 
mance calculations, and in the "Wake” calculation option. The 
variables used for propeller efficiency and wake calculation are 
listed below. 

LOC. 

Default 

Variable 

Description 

1 

0 

WAKBUG 

Debug option. 

1 print WAKEPARM input 
data. 

0 do not print data. 

6 

1 

Y I PLOT 

Wake output. 
1 is on 
0 is off. 

8 

1 

YNPSKN 

Skin friction drag. 


1 : included, 

0 : not included. 


28 


20 

0 

YIPRT 

Output option: 

0 : no output 

1 : V r /V and V x /V 

2 : Fluctuating lift 

23 

0 

YIVWK 

The viscous wake shape is 
described by 1 of 3 op- 
tions: 

0 : gauss pulse 

1 : cosine squared 

2 : cosine 

25 

11 

YNORPE 

No. of output radii. 

751 

.24, .35, YROUTP 
.45, .55, 

.65, .75, 

.80, .85, 

.90, .95, 

.98 

Output radii, r/R^,. 

The following input 

locations may be 

ignored. 

(2) 

.05 

YK 

Spanwise integration mesh 
size, fraction of radius. 

( 9 ) 

.0001 

YTOL 

Fourier series sum con- 
vergence tolerance. 

(10) 

99 

YMM1 

Max. no. of Fourier coef- 
ficients to calculate, 
max. of 99. 

(12) 

.001 

YT0L1 

Tolerance for high fre- 
quency form of Fourier 
coefficient. 

(13) 

0 

YLASTM 

Last Fourier coefficient 
for detailed output. 

(16) 

0.7 

YZO 

Origin of special routine 
for interpolation of cir- 
culation curve. 

(17) 

0.5 

YZNORM 

Circulation curve normal- 
izing factor. 

(18) 

1.0 

YAOPT 

0 to not iterate on in- 
duced angle. 
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( 21 ) 


0 


YMM2 


Turns on diagnostic print 
for Fourier coefficient 
up to YMM2 when YIPRT = 

2 . 
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The following inputs are required in addition to the above if 
wakes are to be calculated. Output radii are divided into those 
within the tip radius, YNOR, and those beyond the tip radius, 
YNPX. 



Description 

Number of output radii 
which are less than or 
equal to 1.0 (within tip 
radius) . 

Number of output radii 
which are greater than 
1.0 (outside of tip). 

No. of output axial loca- 
tions, max of 10 at each 
radius. 

No. of Fourier coeffi- 
cients used in calcula- 
tion of chordwise wake 
component . 

No. of Fourier coeffi- 
cients used in calcula- 
tion of radial wake com- 
ponent . 

No. of Fourier coeffi- 
cients used in calcula- 
tion of downwash wake 
component . 

Output radii, fraction of 
radius. There should be 
YNOR+YNPX of these. 

Output axial locations, 
fraction of radius, for 
1st output radius. There 
should be YXNCS of these. 
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11 


0 


721 0 


YARRY (I,K) 
1=1, YXNCS 


YARRY (I, K) 
1=1, YXNCS 


Output axial locations, 
fraction of radius, for 
second output radius. 
There should be YXNCS of 
these. 

Output axial locations, 
fraction of R^, for third 
output radius. 

There should be YXNCS of 
these. 


For K=1 , YNOR+YNPX 


The following input location may be ignored. 


(54) 0.01 YFPHW Width of averaging func- 

tion for wake calcula- 
tions = X/R^. 


AIRPARMS: 

This command will input options which are used to obtain the 2-D 
airfoil drag from built-in tables of lift and drag coefficients. 

LOC. Default Variable Description 

1 0 AIRBUG Debug Option, 1 to print 


The following input locations may be ignored. 


(2) 24 AIRNUM 


(3) 5 AIRTYP 


(13) 0 AIRCOR 


2-D Airfoil data pack 
number wanted for this 
run. 

Currently, only NACA 16 
is available in this 
code. 

0 : Isolated 2D airfoil 

1 : Isolated 2D airfoil 
w/cascade correction. 

4 : Isolated 2D airfoil 
w/sweep correction. 

5 : Isolated 2D airfoil 
w/sweep and cascade. 

0 : no drag correction 
input 

1 : drag correction 
input 
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III* HUB 



(15) 


0 


CDMLT 


2D drag multiplier, used 
if AIRCOR = 1. 


(17) 0 DCDT 2D drag increment, used 

if AIRCOR = 1. 

CDMLT and DCDT are used to alter the drag of the stored airfoil 
data. They may be used to simulate another airfoil, to account 
for a rough airfoil or provide a better match with test efficien 
cy. The altered drag coefficient is computed as: 

C D(tUered ) = C d *CDMLT + DCDT 
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NOIZPARM ; 

This command will input the options which are used to control the 
acoustic calculations. 


NOTES: B is number of blades, M is noise harmonic order 

B * ZMMAX must be less than 1001 

MX is flight Mach number (may not be 0 for near-field 
calculation) 

LOC. Default Variable Description 


ZNZDBG 


ZNFIND 


ZMMAX 


ALT 


ZNX 


10 


ZXORX1 


Debug option : 

0 no print; 

1 print input parameters; 

2 print file 50 input. 

0 for far-field theory 

1 for near-field theory 

Max noise harmonic order 
group to calculate. 

Distance from prop axis 
to observer, ft. 

Number of axial directiv- 
ities for noise calcula- 
tions, 20 max. 

Axial directivities are 
1 : visual or 
0 : retarded. 


11 XI Visual observer positions 

along axis ( + ahead of 
prop), ft. (calculated if 
X is input) . 


or 
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X 


Retarded positions ob- 
server, ft. (calculated 
if XI is input) . 
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59 1 ZIBLW 


69 90. PHIF 


70 90. PHIR 


The following input locations may 

(51) 10 ZJMINI 

(52) 8 ZNJ1 


1 to add boundary layer 
and wake displacement 
thickness to blade thick- 
ness roonopole noise. 0 
to omit (for comparison 
with other predictions) . 

Azimuthal observer dir- 
ectivity angle for un- 
steady loading, degrees 
to front rotor — normally 
use the default. PHIF 
increases in the direc- 
tion of rotation; zero 
corresponds to the phase 
reference position for 
calculation of unsteady 
air loads. Thus, the ob- 
server position is speci- 
fied relative to the un- 
steady flow field refer- 
ence. 

Azimuthal observer direc- 
tivity angle for unsteady 
loading, degrees to rear 
rotor. Important only 
when BPFF = BPFR. Gives 
the phase difference be- 
tween rotors. PHIR in- 
creases in the direction 
of rotation; zero cor- 
responds to the phase 
reference position for 
calculation of unsteady 
air loads. Thus, the ob- 
server position is speci- 
fied relative to the un- 
steady flow field refer- 
ence. 

be ignored . 

With ZNJ1 , controls 
chordwise integration 
mesh. 

Number of chordwise 
points = ZJMINI + M + B * 
ZNJ1, 1001 max pts. 
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(53) 10 


ZJMIN2 


Used with ZNJ2 , controls 
spanwise integration 
mesh. 


(54) 10 ZNJ2 


(55) 

6 

ZIW 

(56) 

1 

Z START 

(57) 

1 

ZMINC 

(58) 

2 

ZLHSOR 


(60) 

0 

ZPRUNS 

(61) 

0 

ZITRAP 

(62) 

0 

ZNFOUT 


(63) 

0 

ZKK 

(64) 

0 

ZKXTND 

(65) 

10 

ZKMIN 

(66) 

10 

ZOMEGA 

(67) 

0 

ZIPEXT 


Number of spanwise points 
is JMAX (must be less 
than 2001) which is cal- 
culated as follows : 
if ZNJ2 > 0, 

JMAX = ZJMIN2 + 

M*B*ZNJ2 ; 
if ZNJ2 < 0, 

JMAX = ZJMIN2 - ZNJ2 
if ZNJ2 = 0, 

JMAX = ZJMIN2 + ZNJ2 
points/phase cycle. 

Output file number. 

Start harmonic order. 

Increment in harmonic 
order . 

Indicates what type of 
unsteady flow: 

1 = blade wakes; 

2 = flowfield. 

If not 0, print diagnos- 
tic unsteady loading 
noise table. 

If not 0, max diagnostic 
print. 

If not 0, print harmonic 
table to file number 
ZNFOUT; allocate appro- 
priately. 

See note below. 

See note below. 

See note below. 

See note below. 

If not 0, print load di- 
agnostics. 
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( 68 ) 

(71) 

ZKK, 


0 ZIPSIC 


1 ZNAIR 


If not 0, do "instant" 
quadrupole noise calcula- 
tion (not recommended) . 

Code for airfoil thick- 
ness distribution (speci- 
fy for all radial sta- 
tions) : 

1 = series 16 
3 = series 64 
5 = series 65 
7 = biconvex parabolic 
(analytic) 

8=4 digit series 


ZKXTND, ZKMIN, and ZOMEGA : 

In the near-field option, the program evaluates a 
Fourier transform numerically using rectangular 
integration. The range for this frequency integra- 
tion is 1/ ( 1+MX) to 1/ ( 1-MX) , where MX is the flight 
Mach number. The number of steps in the integration 
range determines the trade-off between precision and 
running time. The integrand is typically a fluctuat- 
ing quantity whose rate of oscillation, dl/dw, is 
computed by the program as a function of several 
factors such as harmonic order, Mach number, and 
observer position. To achieve uniform precision over 
a range of conditions, the program determines the 
number of points (ZKK) in the frequency integration 
range from two input numbers: the first is the mini- 
mum points in the range (ZKMIN) and the second is 
the number of mesh points per oscillation of the 
integrand (ZOMEGA) . The number of integration points 
is then ZKK = ZKMIN + ZOMEGA * dl/dw. Default val- 
ues are ZKMIN = 10 and ZOMEGA = 10, which from nu- 
merical tests give a reasonable compromise between 
precision and running time. If the user wishes to 
experiment, he can override either of the two de- 
fault values. 

Reducing the number of points will reduce running 
time and storage requirements, but also reduce the 
precision of the calculation. Determining the satis- 
factory level of tradeoff between precision and 
running time for other than the defaults is up to 
the user. 
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In some cases a significant contribution to the 
noise can be caused by frequencies outside the range 
of frequency integration noted above. To account for 
this, ZKXTND points are included outside both the 
upper and lower bounds of the integration. If ZKXTND 
equals zero then it is automatically computed by the 
program, otherwise the input value of ZKXTND is 
used. However, ZKK + 2* ZKXTND must be less than 401. 
If this is not true then the code will reduce ZMMAX 
by one until this condition is met. 



INTERPRM : 

This command inputs parameters for the control of the non-zero 
harmonic terms in the calculation of the interference between the 
front and rear rotors in counter-rotation cases. These are a 
sub-set of the WAKEPARM parameters. Normally, these default 
values should be used. 

LOC Default Variable Description 

The following locations are optional and may be ignored: 


(1) 

o 

• 

o 

INTBUG 

Printout control: set to 

1 to increase amount of 
printout regarding inter- 
rotor wake. 

(2) 

0.05 

INTXK 

Spanwise integration mesh 
size, fraction of R^. 

(3) 

0.7 

INTZ0 

Origin of special routine 
for interpolation of cir- 
culation curve. 

(4) 

0.5 

INTZNM 

Circulation curve normal- 
ization factor. 

(5) 

1.0 

INTAPT 

0 to not iterate on in- 
duced angle. 

(6) 

1.0 

INTVWK 

The viscous wake shape is 
described by one of three 
options: 

0: Gaussian pulse 

1: cosine squared pulse 

2: cosine pulse 

(7) 

0.01 

INTFPH 

Width of averaging func- 
tion for wake calcula- 
tion, in terms of R,^. 

(8) 

1.0 

INTSKN 

Skin friction drag. 
1: included 

2: not included 

(9) 

0.001 

INTTOL 

Fourier series sum con- 
vergence tolerance. 

(10) 

0.0001 

INTOL1 

Tolerance for high fre- 
quency form of Fourier 
series . 
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(11) 

99.0 

INTMM1 

Maximum number of Fourier 
coefficients to calcu- 
late, limit of 99. 

(12) 

0.0 

INTMMU 

Number of Fourier coeffi- 
cients used in calculat- 
ing the chordwise compo- 
nent of the wake. 

(13) 

0.0 

INTMMV 

Number of Fourier coeffi- 
cients to be used in cal- 
culating the radial por- 
tion of the wake. 

(14) 

0.0 

INTMMW 

Number of Fourier coeffi- 
cients to be used in cal- 
culating the downwash 
portion of the wake. 

(15) 

0.0 

INTPLT 

Wake output; 1 is on; 0 
is off. 

(16) 

0.0 

INTPRT 

Output option. 


0 : no output 

1: print V r /V and V x /V 

2: print fluctuating 

lift. 
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VORTPARM : 

This command inputs the options which are available to control 
the Vortex Flow Aerodynamic calculations. This command has two 
sub-commands: FWD and AFT. Use VORTPARM (FWD ) to enter 

parameters for the front rotor in a counter-rotation case. Use 
VORTPARM (AFT )* to enter parameters for the rear rotor in 

counter-rotation cases. Use VORTPARM, with no sub-commands, for 
single rotation cases. 


LOC. 

1 


2 


3 


4 


5 

6 


Default Variable 

0 VTXDBG 


1 DOLEAD 


1 DOSIDE 


1 DOAUGL 


.97 ZAUGFF 

0 TIPLOR 


Description 

Debug option: 
if » 0, input not 
printed, 

if = 1, input printed. 

If = 1, then calculate 
additional lift due to 
leading edge vortex. 

If = 1, then calculate 
additional lift or 
radial force due to tip 
edge flow. 

If = 1, then calculate 
additional lift due to 
leading edge vortex 
shed over aft portion of 
blade at the tip 
(augmented lift) . 

Radius at which augmented 
lift acts, r/R^. 

Indicates whether tip 
edge flow results in ex- 
tra lift (0) or radial 
force (1) . 


TIPLOR is set by the user to determine the type of 
edge flow. 

0 for separated tip flow, where the tip vortex 
gives extra flow at the tip. 

1 for attached tip flow which produces a radial tip 
edge force. 


*The AFT and FWD subcommands are eight characters long. The 
trailing blanks are not required. 
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The following input locations may be ignored: 


(7) 0 PRTPRF 


(8-57) 0.2, VTXDLX 

at all 
radii . 


1: print results on ev- 

ery iteration. 

0: do not print results 

on every iteration 

Chord-wise distance from 
leading edge at which 
vortex loads are assumed 
to act, in terms of frac- 
tion of local chord. 

There is one value of 
VTXDLX for each output 
station in the spanwise 
direction, specified as 
ZNIS (location 950) in 
the 2-D coordinate inputs 
under blade geometry 
(BLADEGEO section above) . 
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AEROEXEC ( PRNTCASE ) : 

This command (sub-command) will process all of the input and 
execute the code to the point where the Compressible Panel Method 
prints the data it has received. The aerodynamic calculations 
will not be performed. This is useful in verifying the input 
before doing the actual aero calculations. 

AEROEXEC (EXECCASE) ; 

This command (sub-command) will process all of the input and 
execute the aerodynamic calculations. 

NOIZEXEC (EXECCASE) ; 

This command (sub-command) will execute the noise calculation 
section of the code. This should be placed after the AEROEXEC 
record . 

WAKEEXEC (EXECCASE) : 

This command (sub-command) will execute the potential and viscous 
wake calculations. This should be placed after the AEROEXEC or 
NOIZEXEC record. This command should only be used with single 
rotation operation of the code. 

ENDCASE ; 

This command signifies the end of all commands for the current 
case, and is required. 

END JOB : 

This command will terminate execution, and is required. 
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V Output Description 

This section presents a description of the output from the UCAP 
program and includes output print from AEROEXEC (the panel code 
performance calculation) . The output from NOISEXEC (the noise 
calculation) is described in Section VI. 

The output of this program expands on the output of the UAAP 
program (reference 3) . The outputs of the blade geometry rou- 
tines, the vortex routines, and the lifting surface (steady) 
solution routines have changed due to the improvements in the 
code: these changes are applicable to single and counter- 

-rotation operation. Also, the program output has been expanded 
due to the addition of counter-rotation; the outputs for stream- 
line contraction, rotor-to-rotor interference velocities, and 
rotor control point location are new. 

The debug option, described under the various input section 
commands, controls the scope of output. These options normally 
have values of either 0.0 (the default), for minimal printout, or 
1.0, yielding additional output. Most pages have a title fol- 
lowed by the header information (entered via the HEADER command) . 
The header consists of a brief description of the page contents. 
This includes, after the first colon, the subroutine name which 
is doing the printing of this page, 2) the input command option 
currently in effect, 3) the time and date at the start of the 
run, and 4) the program/version identification. 

Output Description when DEBUG = 1. 

This output description follows the order shown in Section IV, 
"Input Description". 

OUTPUT ; 

The output starts with an echo of the input data set, which is 
discussed in the input description, and shown in Figure 11. 

RUNPARMS : 

The debug option of 1. will print the page shown in Figure 12, 
which shows the input parameters, defined in the input section, 
along with the location number and the input or defaulted value. 

CRPPARMS : 

The debug option of 1 will print the page shown in Figure 13, 
which shows the parameters, defined in the input section, to be 
used for counter-rotation. FWD->AFT HARM, and AFT->FWD HARM 
specify the number of unsteady interference harmonics for the 
front rotor upon the aft rotor and the aft rotor upon the front 
rotor, respectively. 
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AIRPARMS : 

Figure 14 presents the additional output for debug = 1. Again, 
the parameters are defined in the input section. A few lines of 
airfoil description are given, defining the airfoil selected, 
used to determine profile drag. 

BLADEGEO : 

Setting DEBUG =1. in the BLADEGEO section of input yields an 
additional page of input definition, given in Figure 15. The 
sample shown in Figure 15 is for the forward rotor of a counter- 
rotation Prop-Fan. The items output and the form of the output 
is the same for either rotor of a counter-rotation Prop-Fan or 
the rotor of a single rotation Prop-Fan, so samples for those 
cases will not be included. Definitions of the parameters 
follow: 


PARAMETER 

DEBUG 

BLADN 

DIAMETER 

SCO 

SWEEP TYPE 
DESIGN ANGLE 
RUNNING ANGLE 
TYPE CUT 
X 

T/B 

B/D 

CAMBER 
DELTA THETA 

XSWP 

YSWP 

ZSWP 


DESCRIPTION 

Debug option for additional printout. 
Number of blades. 

Propeller diameter, feet. 

Inner-most blade station at 50% chord. 
Defined in input as SWPOPT. 

Defined in input as THTDES . 

Defined in input as THTCUT . 

Defined in input as ZKCUT. 

Station radius/blade tip radius (R^ p ) . 
Blade thickness/blade chord. 

Blade chord/propeller diameter. 

Design lift coefficient. 

Blade twist, from plane of rotation, de- 
grees. 

Defined in input. 

Defined in input. 

Defined in input. 
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AIRFOIL TYPE 

Defined 

in 

input 

as 

BAFL. 

NO. OF % CHORD 

Defined 

in 

input 

as 

ZNPCOV . 

LIST OF % CHORD 

Defined 

in 

input 

as 

PCTCHD . 

NO. OUT STATIONS 

Defined 

in 

input 

as 

ZNIS. 

LIST - STATIONS 

Defined 

in 

input 

as 

ZBLDST. 


LSTPARMS ; 

Selecting the debug option, LSTBUG = 1, will print additional 
input definition as shown in Figure 16, which lists the available 
input options, their location number and the selected or default- 
ed values. Refer to the input description for more detail on 
input items. This is repeated for both rotors of a counter- 
rotation Prop-Fan. 


NOIZPARM : 

The debug option of 1 will print the page shown in Figure 17, 
which shows the input parameters, defined in the input section. 


VELGRADS : 

The debug option of 1 will result in the printout shown in Figure 
18. Note that this input is not from the same test case as most 
of the other output sections and the positions and velocities are 
not normalized. 


VORTPARM : 

The debug option of 1 will print the page shown in Figure 19, 
which shows the input parameters, defined in the input section. 
The form of the VORTPARM output is identical for both the front 
and rear rotor in counter-rotation cases. Values of "X FOR LEV 
ACTION" are the chordwise positions (VTXDLX) at which the vortex 
loads act (see inputs) . 


WAKEPARM : 

For counter-rotation operation, no items should be entered under 
the WAKEPARMS command: this description is superfluous. For 

single rotation, the printout in Figure 20 will be obtained. The 
input parameters are defined in the input section. 
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TNTERPRM ; ^ , * 

If the debug parameter for INTERPRM is set to 1, the output in 

21 results. The input parameters are defined m the input 
section. The following table relates the input and output 
designations. 


OUTPUT ITEM 

CORRESPONDING INPUT ITEM 

Debug Switch 

INTBUG 

K 

INTXK 

ZO 

INTZO 

ZNORM 

INTZNM 

AOPT 

INTAPT 

IVWK 

INTVWK 

FPHW 

INTFPH 

NPSKN 

INTSKN 

TOL 

INTTOL 

TOLI 

INTOLI 

MMI 

INTMMI 

IPLOT 

INTPLT 

IPRT 

INTPRT 
























Aerodynamic Output 


AEROEXEC i 

Figures 22 through 28 show the printout when LSTPARMS and 
CRPPARMS debug variables are set to one. Figure 22 shows the 
output for the GETXMB routine, which calculates the locations of 
the control points on each rotor relative to the other. This 
routine prints out seven blocks of information: some general 

propeller information, the geometry information for the forward 
rotor, the location of the forward rotor control points relative 
to the rear rotor, the forward rotor leading edge sweep informa- 
tion, the geometry for the rear rotor, the location of the rear 
rotor control points relative to the forward rotor, and the rear 
rotor leading edge sweep information. 

Figure 23 shows the blade description. MCA and FA are shown in 
Figure 24. These describe the location, in space, of the mid 
chord line of the blade. 

Figure 25 shows the locations of the front row control points 
relative to the pitch change axis of the rear rotor. DELTA PHI 
is the angle measured from the helicoidal surface to the control 
point. XMBAR is the distance, parallel to the axis of rotation, 
from the pitch change axis to the control point. MBAR is the 
control point number at a given radius. MU is the control point 
number counting from the forward control point at the root to the 
aft-most at the tip. 

Figure 26 is the list of forward rotor inner and outer wake 
points. The inner wake points are those within a circular 
cylinder defined by the tips of the rear rotor. Outer wake 
points would be outside this cylinder. 

Figure 27 is output of the LESWP routine, which calculates 
leading edge sweep. The sweep angle is shown in degrees at the 
stations specified by ZNIS and ZBLDST in the BLADEGEO inputs. 

Figure 28 is the list of rear rotor inner and outer wake points. 
The inner wake points are those within a circular cylinder 
defined by the tips of the forward rotor. Outer wake points 
would be outside this cylinder. 

Figure 29 shows the re-mapping of each rotor's control points, 
relative to the other rotor assuming the upstream influence 
follows streamlines. Since no data were entered under VELGRADS 
in this case, the default number of 5 streamlines (4 streamtubes) 
is used, the streamlines are concentric cylinders and the re- 
mapping does not move anything. Note that the first control 
radius is below the spinner cut-off, and therefore outside the 
inner streamtube (thus has the null streamtube 0 assigned to it) , 
and that the outer 4 control radii are all contained in the outer 
streamtube (streamtube 4) . 
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Figure 30 shows 
solver (F271) . 

PARAMETER 

TEMP, DEGF 

RHO/RHO STD 

SPEED OF SOUND 

ADVANCE RATIO 

FLIGHT MACH NO. 
FLIGHT SPD KTS 

RPM 

TIP HEL. MACH 

START BLENDING 
DIAMETER 
NO. BLADES 
NO. INPT STA. 
FREQ. OF UNST. 
NO. NODAL DIA. 

K - DOWN 
K - START 
INPUT STATIONS 
B/D 

TOTAL TWIST 
MCA/D 


the information passed to the lifting surface 
This is: 


DESCRIPTION 

Ambient temperature, °F 

Ambient density/ standard density 

Ambient speed of sound, fps 

V/nD, where V is flight velocity, in fps, 
n is propeller rotation rate in revolu- 
tions/second, and D is diameter, in feet. 

Flight Mach number 

Flight velocity, knots (nautical miles per 
hour) 

Propeller rate of rotation, rpm 

Mach number based on resultant of flight 
and tip speeds. 

Defined in LSTPARMS input as QMBLEN 

Propeller diameter, feet 

Number of blades 

Number of input stations 

Defined in LSTPARMS input as QQ 

Defined in LSTPARMS input as QK 

Defined in LSTPARMS input as QKDOWN 

Defined in LSTPARMS input as QK START 

Blade radial station/blade tip radius 

Blade chord/ propeller diameter 

Operating twist from plane of rotation, 
degrees 

Distance from pitch change axis to blade 
mid-chord point along the helix/propeller 
diameter, see Figure 24. This and FA/D 


49 





(next item) are functions of blade geome- 
try and advance ratio. 

FA/D 


Perpendicular distance from helix to the 
intersection of the blade mid-chord and 
mid-camber point/diameter, see Figure 24. 

T/B 


Maximum blade thickness/blade chord 

CLD 


Blade design lift coefficient 

SWEEP 


Sweep angle between mid-chord line and 
resultant inflow velocity 

LSI CAMPER TABLE: 

Another output table is shown in Figure 31. This is a table of 
the blade camber angle, measured from the plane of rotation as a 
function of local blade radius (normalized by tip radius) and 
fraction of the local chord. 


ANSO : ■ 

Figure 32 gives the initial output from subroutine ANSO. This 
includes the contributors to the velocity triangle at each radial 
station and each control point. There are eleven columns; they h 

are describe below. m 


Parameter 


Description 


I,MBAR, MU 


R/R 

TH+CAMB 


Indexing information. I is 
station number, counting from 
root; MBAR is control point 
number, counting from leading 
edge of blade. MU is control 
point index, counting from lea- 
ding edge to trailing edge, 
from root to tip. 

Radial location of station in 
terms of R/R,^. 

Angle from plane of rotation to 
camber surface of the blade at 
this location, measured in 
degrees. 
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THICKNESS 


V-AXIAL 


V-SWIRL 


PI*X/J 


CENTER BODY 


WMU 


W T , thickness vector or the 
contribution to turning due to 
airfoil thickness, measured in 
degrees. In this case the W T 
could not be calculated due to 
numerical problems on the IBM 
3090. Typically the magnitude 
of this is less than 2°. 

Axial component of the local 
change in velocity, AV/V 0 , due 
to induced flow and the pres- 
ence of another rotor. This 
will be zero for the first pass 
for the forward rotor. 

Tangential component of the 
local change in velocity, 

AV/V 0 , due to induced flow and 
the presence of another rotor . 
This will be zero for the first 
pass through the forward rotor. 

7T • ( r /Riip) / J , the tangent of the 
advance angle at each radial 
station. 

The local change in velocity 
due to the presence of the cen- 
ter body; this is from the 
velocity field which has been 
entered via the VELGRADS com- 
mand. This is zero as no ve- 
locities were entered via the 
VELGRADS command. 

Turning required to ensure the 
flow is tangent to the camber 
surface. 


Also printed here is the coefficient of sound power, which is; 
acoustic power loss/ (density* (rev/unit time) 3 * (diam) 5 ) 


PERFORMANCE RESULTS : 

Figure 33 shows the performance results for the individual rotor; 


C p =P/pn 3 D 5 
C t =T/pn 2 D 4 

Efficiency = J(C,/C p ) 
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where 


n is rotor rotational speed 
D is diameter 
p is fluid density 
P is power 
T is thrust 

The units of n, D, and p, P, and T must be consistent. 

Included are potential loads plus profile drag losses and the 
various non-linear lift terms which may result from leading edge 
sweep . 

3&RTEX CALC. ; 

Figures 34 through 38 show the vortex calculation outputs. 

Figure 34 shows the operating conditions. Figure 35 shows the 
angle of attack (A.O.A) at the leading edge, three dimensional 
angle, (ALPHA 3-D) , induced angle, advance angle, and blade 
angle. Figure 36 shows the coefficients from the potential 
calculations. CDPOT does not include leading edge thrust. L.E. 
K, MAG is the leading edge suction force coefficient. This is 
calculated within the lifting surface portion of the code. It is 
used to calculate the leading edge force coefficient. Figure 37 
shows the lift and drag coefficients which result from the vortex 
calculations. The incremental lift and drag coefficients illus- 
trate the relative magnitude of the vortex components. Note that 
the tip edge vortex can cause lift (separated flow) or radial 
forces (attached flow) . The terms "TIP VORTEX" and "SE VORTEX" 
(side edge) are used interchangeably. Figure 38 shows the 
elemental performance. MODOPT and NBOPT are LSTPARMS inputs and 
are explained in the LSTPARMS input section. 

PERFORMANCE RESULTS : 

Figure 39 shows the performance of each individual rotor. CP is 
the standard propeller power coefficient, CP=P/ (pn 3 D s ) . CT is the 
standard propeller thrust coefficient, CT=T/ (pn 2 D 4 ) . This infor- 
mation is repeated for each pass through each rotor. 


I NTERS : 

Figure 40 shows the velocities calculated as a result of the 
steady interference between rotors. The first section is a 
tabulation of the change in axial velocity due to the center body 
(VX-BODY) , which is zero for this case as no velocities were 
entered with the VELGRADS command, and the rotor's induced 
velocity. VX-POT is the change in axial velocity, AV X /V 0 due to 
potential flow. VX-2DARY is the increment due to non-potential 
loading, e.g., vortex flow. The second section shows the veloc- 
ity increment due to rotor loading at the rear rotor control 
points. 


52 


ksA 




PRTPRF • 

Figure* 41 shows the performance summary at one of the intermedi- 
ate steps. In part A of this Figure the performance without the 
vortex loading is shown, in part B the performance with the 
vortex loading. Both parts are divided into two sections: a 

coefficient section and a ratio section. In the coefficient 
section, the power, thrust, and torque coefficient and the ef- 
ficiency are shown. The last coefficient, ROLLING MOMENT, is 
defined as 

C rm = (Qfwd“Q»ft) (pn 2 D 5 ) 


where Q^., and Q-ft are the forward and rear rotor torque, respec- 
tively. All coefficients in the "TOTAL" column are based on the 
diameter and rotation rate of the forward rotor, hence FWD + AFT 
do not equal TOTAL. The ratio section has the ratio of tip 
speed, diameter, and angular speed of the aft rotor to the 
forward rotor. Also given are ratios of area, power, thrust, and 
torque for forward rotor to total, aft rotor to total, and aft 
rotor to forward rotor. 


Figure 42 shows the output of the routine which performs stream- 
line adjustment. The first section shows the axial velocities at 
each control point. The second section shows the flow rates 
through each stream tube. The third section ^ shows the new 
boundaries to the stream tubes after continuity is enforced. 


Figure 43 shows the re— location of the control points which are 
used for calculation of the flow induced by the rear rotor on the 
forward rotor. 


CPMAP: 

Figure 44 shows the relocation of the control points which are 
used for the calculation of the wake flow from the front rotor 
upon the rear rotor. 

The output produced for each iteration is not included here; 
there is about 26 pages of printout per iteration with the DEBUG 
variable set to 1 in all input sections. Due to volume, they are 
not included here. 

PROPELLER EFFICIENCY PROGRAM : 

Figure 45 and 46 present the rotor-to-rotor interference field 
calculated by the near wake module. The information in Figure 45 
is the steady calculation. At each station, there are two sets 
of data. The first is a single line with: 

STATION Station number, (1 = innermost) radial 

location in terms of R,^ of radius gener- 
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PHI 


SIGMA1 

SIGMA2 

SIGMA1/SIGMA2 

TWIST 


ating influence field (here, the forward 
rotor) 

Advance angle, degrees 

Internally calculated values used in 
calculations of wake velocities which 
depend on radius and operating condition 

Angle between chordline and plane of rota- 
tion 


Next there are ten lines (One line for each of QNCP points at 
each radius) 


MBAR 

XMBAR 

WMU 

(STEADY) 

DOWN WASH 

CHORDWISE 

VISCOUS 

WAKE 


Chord-wise control point number 

Axial distance from pitch change axis, 
normalized by tip radius of rotor generat 
ing influence field 

Perturbation angle produced by field 


Down wash component of WMU 

Chord-wise component of WMU 

Portion of WMU produced by the viscous 
deficit caused by profile drag of the 
front rotor 


AXIAL 


Axial component of WMU 


SWIRL 


Swirl component of WMU 


Figure 46 shows a sample printout for the higher harmonic portion 
of the rotor-to-rotor interference field. 

The first five lines show: 

Direction of interference 
Harmonic of blade pass frequency, m 
Number of blades, B 

Harmonic of shaft frequency, described in Section IV, "Input 
Description", under LSTPARMS as QK, and calculated as: 

K f r oot =m * B »ft 

K ift =m *Bf r0 o t 
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Q-order, described in Section IV, "Input Description", under 
LSTPARMS as QQ, and calculated as: 

Qfrool = ^froot~ (i"*" ( *"*froot/ ^»ft) ) 

Qaft =K »ft“ ( (^aft/^froot) ) 

K and Q are calculated for the rotor receiving the interference. 

Next, is one line for information for items which are constant at 
a given radius: 

STATION Station number and distance from center of 

rotation normalized by R^. 

WAKEPNT This is the normalized contracted radius 

for this station. 


PHI 


Advance Angle, degrees. 


SIGMA1 , 

SIGMA2 , 
SIGMA1/ SIGMA2 


Internally calculated values used for cal- 
culations in the wake. These depend on 
radius and operating condition. 


TWIST 


Angle between chordline and plane of rota 
tion, degrees. 


After the line for each radius, there is one line for each 
chordwise point (QNCP) : 

MBAR, MU Indexing information. 

XMBAR Distance from pitch change axis of rotor 

generating interference field. 


DPHI Angle from helicoidal surface to control 

point, degrees. 

WMU Unsteady loading downwash vector. This is 

the QWMU value described under LSTPARMS in 
Section IV, "Input Description". 


DOWNWASH 


Component of WMU normal to camber line of 
blade at this radial station. 


CHORDWISE 


Component of WMU tangent to chordline at 
this radial station. 


VISCOUS Component of WMU due to viscous wake. 

WAKE 
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VI. Noise Outputs 


NOZOUT: 

The first page of the noise program output (see Figure 47) 
summarizes the observer positions selected by the user. The 
first few lines repeat the header information. Then, the computer 
time used is presented. Then, the observer locations selected by 
the user are printed. First, the sideline distance for the front 
(or single) and/or rear rotors is printed. Then, for the front 
rotor, the retarded axial positions, the visual axial positions, 
and the corresponding retarded and visual radiation angles are 
printed. The relationships between the visual and retarded 
distances and angles are also shown in Figure 47. Finally, the 
front rotor azimuthal observer angle is printed. 

The rear rotor is offset axially from the front rotor by a 
fraction of the front rotor radius. This fraction is printed, and 
the corresponding retarded and visual axial positions and radia- 
tion angles are printed. Listing of the rear rotor azimuthal 
observer position completes the observer geometry printout in 
figure 47. 

The local ambient conditions are then listed, giving the tempera- 
ture, ambient pressure normalized by sea level standard, and 
local speed of sound. 

On the next page, in Figure 48, the front and rear rotor operat- 
ing conditions are printed. FLIGHT MACH NUMBER is the free- 
stream Mach number. TIP ROTATIONAL MACH is the tip speed 
divided by the speed of sound. TIP RELATIVE MACH NUMBER is the 
resultant of the flight and tip rotational Mach numbers. MT/MX 
is the ratio of tip Mach number to flight Mach number. BPF is 
the blade passing frequency, measured in Hertz. The 
FOOT DIAMETER is tip diameter, in feet, of the rotor. 

On the next page (Figure 49) is a header indicating the beginning 
of the noise calculation results. Following the header is a 
summary of the noise calculation results for the first directiv- 
ity. This type of page is repeated for as many directivities as 
were requested by the user. The X POSITION ANGLE for visual and 
retarded axial position and angles for front and rear rotor is 
shown. Then, for each possible noise radiation frequency, a 
summary of the noise radiated by the front and rear rotors is 
given. The M COUNTER and K COUNTER are given for the front rotor. 
Thus, K COUNTER = 0 corresponds to steady loading, and the steady 
sources are also printed. The roles of the M COUNTER and 
K COUNTER are reversed for the rear rotor, so M COUNTER = 0 
corresponds to steady loading for the rear rotor. Results are 
presented for each harmonic up to MMAX, and then a new page (not 
shown) presents results for the next axial observer location. 
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Finally, a summary of the total noise at each frequency is pre 
sented in Figure 50. The header informatiort is followed by a 
listing of the observer coordinates. The computed radiation 
frequency and the total combined noise for the front and rear 
rotors are presented at each frequency. This output format is 
repeated for all requested directivities and harmonics. 


VII. Program Installation 

UCAP is written in FORTRAN-77, mostly adhering to ANSI X3. 9-1978, 
with the common additions of the double precision complex 
(complex*16) data type and the IMAG intrinsic function. UCAP 
requires the IMSL Libraries, version 1.1 (formerly IMSL 10), 
mostly for FFT's, Bessel functions, and matrix inversions. The 
IMSL Libraries are copyrighted; they are not supplied with UCAP. 

A list of the routines directly called from UCAP is provided in 
Table II. 

The target operating environments for this program are the Cray 
series of super-computers. This program was developed on an IBM 
3090 running MVS/XA and a Sun workstation running Sun's version 
of UNIX. Based on experience with the interim (Task 2) version 
of UCAP and with UAAP, this program is expected to be operated on 
the Cray with double precision turned off (by a compiler option) 
and with the single precision version of the IMSL routines. Many 
complex numbers are operated on within UCAP; Table III presents 
a list of subprograms where AIMAG or IMAG is used. In the IBM and 
Sun dialects of FORTRAN-77, the IMAG intrinsic function is 
required when operating on complex*16 data. This is replaced by 
AIMAG for Cray FORTRAN, which does not support the complex*16 
data type. Because of the different IMSL calls and the need to 
use AIMAG, directives for cpp (the C pre-processor, a utility 
supported by UNICOS and most UNIX systems) were included into the 
distribution source code to select single precision IMSL calls 
and AIMAG on Cray systems versus double precision IMSL calls and 
IMAG on other systems. 

On Cray systems running UNICOS, selection of single/double 
precision IMSL subroutine calls and I MAG/ AIMAG intrinsic function 
use would be done by compiling with: 

cf77 -DCRAY -c -dp ucap.F 

if the FORTRAN compiler supports the C pre-processor or 

cpp -DCRAY ucap.F > ucap.f 
cf77 -dp -c ucap.f 

where the FORTRAN compiler does not. 

On Cray systems, segldr would then be used to link the resulting 
object code (ucap.o) with the IMSL Libraries in order to produce 
an executable copy of the program. 
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VIII. Pile Requirements 

UCAP requires several files for input, output, and data transfer. 
The input data file (number 5) , and the input data files for the 
noise routines (numbers 49 and 50) are fixed record length files. 
These must have a minimum of 80 characters per record, excluding 
the new line or carriage return used to terminate lines in the 
UNIX system; these will have to be padded with blanks. The 
input data files for the noise routines generated by UCAP are 
properly formed for program operation. 


File Number File Contents 


File Remarks 


5 


6 


Input Data 


Printed Output 


Fixed record length, with a 
minimum of 80 characters per 
record. A typical input file 
will contain about 400 lines 
(32,000 bytes) . 

Standard FORTRAN printed out- 
put: up to 133 characters per 

record, with the first charac- 
ter reserved for carriage 
control. This will be about 
10000 lines. 


11 


40 


41 


The following files are generated by the program: 


Copy of Input 
Data 


Internal Storage 


Internal Storage 


Fixed record length, 80 char- 
acters per record. This is a 
"scratch" file which need not 
be retained after program 
execution. This will be the 
same size as file used for 
unit number 5. 

"Scratch" file for internal 
use. This must be about eight 
million bytes (one million 
Cray words) . 

"Scratch" file for internal 
use. This must be about eight 
million bytes (one million 
Cray words) . 
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49 


50 


51 


52 


AFTTHV 

FWDTHV 

SRPTHV 


Aft Rotor Input 
Data for Noise 
Routines 


Forward Rotor 
Input Data for 
Noise Routines. 


Internal Storage 


Internal Storage 


Aft Rotor Thick- 
ness Vector 
Forward Rotor 
Thickness Vector 
Single Rotor 
Thickness Vector 


Fixed record length, 80 char- 
acters per record. This will 
be about 400 lines for the 
steady and 200 additional 
lines for each harmonic speci- 
fied in CRPPARMS. 

Fixed record length, 80 char- 
acters per record. This will 
be about 400 lines for the 
steady and 200 additional 
lines for each harmonic speci- 
fied in CRPPARMS. 

"Scratch” file for internal 
use. This must be about eight 
million bytes (one million 
Cray words) . 

"Scratch" file for internal 
use. This must be about eight 
million bytes (one million 
Cray words) . 

Unformatted, i.e., written 
without FORTRAN format con- 
trol. These may be saved for 
later use. Size is about 800 
bytes (100 Cray words) . 


Due to the necessity of generating numerous influence coefficient 
matrices, UCAP uses the FORTRAN open statement to assign file 
names at run time, when running counter-rotation. Any influence 
coefficient matrix file may be retained for use in later program 
runs. However for counter-rotation the program must re-generate 
all influence coefficient matrices for a given rotor if any new 
influence coefficient matrix is required for that rotor. 


File Name 


File Contents 


File Remarks 


SRP000 


AFT000 

FWD000 


Single rotation influence 
coefficient matrix. This 
is used for steady and 
unsteady matrices. 


Unf orma tted . Ap- 
proximately 20,000 
Cray words (160,000 
bytes) . 


Counter-rotation steady 
influence coefficient 
matrices, for the aft or 
forward rotor. 


Unformatted. Ap- 
proximately 20,000 
Cray words (160,000 
bytes) each. 
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AFT001 through 
AFTOnn 


Counter-rotation unsteady 
influence coefficient 
matrices, for the aft 
rotor. The trailing 
number (Onn) is incre- 
mented by one for each 
unsteady matrix required 
for the number of inter- 
action orders (forward to 
rear rotor) specified in 
CRPPARMS . 


Unformatted. Each 
is approximately 
20,000 Cray words 
(160,000 bytes) 
each. 


FWD001 through 
FWDOnn 


Counter-rotation unsteady 
influence coefficient 
matrices, for the forward 
rotor. The trailing 
number (Onn) is incre- 
mented by one for each 
unsteady matrix required 
for the number of inter- 
action orders (forward to 
rear rotor) specified in 
CRPPARMS. 


Unformatted. Each 
is approximately 
20,000 Cray words 
(160,000 bytes) 
each. 
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IX. Error Messages 


This section contains a description of the printed error 
messages, along with possible causes and remedies in alphabetized 
order . 

"AIRFOIL NO. = ... OFF AIRFOIL DATA AND IOFF ..." 

The program obtains profile drag from airfoil data 
tables as a function of lift, mach number, airfoil 
thickness ratio, and airfoil chord. This message 
indicates one or more of the above parameters ex- 
ceeded the limits of the stored data. 

Recommendation : check the output for reasonable 
values of drag coefficients. Use AIRPARMS inputs 
CDMULT and DCD to correct unreasonable values. 

"AIRFOIL . . . OFF AIRFOIL DATA " 

Same as above. 

"BDSO ..." 

There are a number of errors, generally resulting in 
program termination, which are associated with the 
blade geometry generator; these messages are pref- 
aced with BDSO. They arise because the blade surface 
generated from user input in BLADEGEO was not suffi- 
cient to allow intersection with the output stations 
(conical surfaces) specified in BLADEGEO variable 
ZBLDST. 

Solution: The problem can usually be avoided (with 
only a slight loss in program accuracy) by setting 
THTDES = THTCUT in BLADEGEO, and assuring that X(l) 

< ZBLDST (1) and X(10) > ZBLDST (ZNIS) in BLADEGEO 
input. 

••*** CONTROL POINT STATIONS NOT WITHIN ..." 

The radial difference between control points, QZAR 
in LSTPARMS, and ZBLDST in BLADEGEO > .02 . 

Solution : Change the either QZAR or ZBLDST input 
values. 

"FAILED TO CONVERGE ON CL ..." 

An iteration failed when trying to obtain the pro- 
file drag data from the airfoil tables. 

Recommendation : check the output for reasonable 
values of drag coefficients. Use AIRPARMS inputs 
CDMULT and DCD to correct unreasonable values. 
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"FOR OUTPUT POINT NUMBER ..." 

Interpolation of the flowfield from VELGRADS is 
performed in two ways : radial first, axial second, 
and then axial first, radial second. Usually these 
interpolations yield the same value of velocity at a 
fixed radial and axial location. However, for this 
point the difference in interpolations exceeded 1 %. 
Probable cause : the flowfield is not smooth in 
either or both the radial and axial directions. 

"FOR OUTPUT POINT NUMBER ... THE AXIAL ..." 

The axial extent of coordinates input into VELGRADS 
was not sufficient to allow interpolation of the 
flowfield at either the blade leading or trailing 
edge. 

Solution: add more points to the flowfield input in 
VELGRADS in the axial direction. 

"FOR OUTPUT POINT NUMBER . . . THE RADIAL ..." 

The radial extent of coordinates input into VELGRADS 
was not sufficient to allow interpolation of the 
flowfield at either the blade root or tip. 

Solution : add more points to the flowfield input in 
VELGRADS in the radial direction. 

"INVALID INPUT CHARACTER ..." 

The load routine has found a character in column 1 
which is not an "L","C","E", or blank. 

Possible causes : "END" record omitted; numeric data 
in column 1. 

"THE INPUT CONTAINS AN UNKNOWN COMMAND ... " 

A COMMAND was expected at this point in the input, 
however, the 8 printed characters do not represent a 
recognizable command. 

Possible causes : misspelled command; no "END" re- 
cord in load format, incorrect case (use UPPER CASE 
only); COMMAND didn't start in column 1. 

"THE INPUT CONTAINS AN UNKNOWN SUBCOMMAND.." 

A sub-command was found which didn't match the ac- 
ceptable sub-commands. 

Possible causes : incorrect spelling; incorrect 
placement of parenthesis. 

"THE NUMBER OF HEADER CARDS " 

More than 10 header records were encountered during 
processing of the HEADER command. 

Possible causes : more than 10 records following the 
HEADER command; no "END" record after the last head- 
er record. 
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"V1S1X1 DIMENSIONED 20X20 ..." 

The flowfield input in VELGRADS was too large; the 
maximum number of coordinates in the axial or radial 
direction is 20. 

Solution : reduce the number VINRD and/or VINAS to 

20 . 

****** MM2 AND MM3 MUST BE LESS THAN 100 ...» 

See QMM2 and QMM3 limits in LSTPARMS . 

"***** NBOPT MUST BE 1 OR 2 ...» 

See QNBOPT limits in LSTPARMS. 

"***** NCP GREATER THAN MAX NCP ...» 

See QNCN limits in LSTPARMS. 

"***** NIS GREATER THAN MAX NIS ... » 

See ZNIS limits in BLADEGEO. 

»** NOIZEXEC CALLED WITH KEYWORD . . T » 

An unrecognizable sub-command was found. The sub- 
command must be "EXECCASE". 

»•***** NSM GREATER THAN MAX NSM ... » 

See QNSM limits in LSTPARMS. 

••**«** NSM * NCP GREATER THAN ...» 

Ensure that QNSM * QNCP < 1000 in LSTPARMS. 

»*** NUMBER OF NODAL DIAMETERS (K) MUST ...» 

If QQ = 0 then QK must equal zero in the LSTPARMS 
input. 

****** N x GREATER THAN ...» 

See ZNPCOV limits in BLADEGEO. 

••**** THE DATA FOR THIS RUN & THE DATA ...» 

A label is attached to the "K-INVERSE" matrix de- 
scribing the parameters used to create it. The cur- 
rent value of the input variable listed below the 
error message does not agree with that in the at- 
tached label. 

Cause : This generally arises because a "K-INVERSE" 
matrix from a previous run was used as input to this 
run (see Q PARTI in LSTPARMS) and the data used to 
generate that matrix is different than that being 
used to run the current case. 

••*** THE NO. OF DIRECTIVITY POINTS ...» 

See ZNX limits in NOIZPARM. 
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***** THE no. of harmonics is greater ..." 

See ZMMAX limits in NOIZPARM. 

***** THE NO. OF HARMONICS IS GREATER THAN ALLOWED 150" 

***** THE NO. OF HARMONICS X THE NO. OF BLADES IS LIMITED TO 
1000 " 

***** THE NO. OF RADIAL INTEGRATION POINTS IS LIMITED TO 2000" 
***** THE NO. OF DIRECTIVITY POINTS IS LIMITED TO 20" 


Cause - input values larger than permissible 
Solution - correct input 
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Appendix A 

Example QWMU Generating Program and Modified Block Data CRPBLK 

This section has the listing (Figure 51) for a small, stand alone 
program, QWMUGEN, which will generate the unsteady downwash 
vector, QWMU, described under LSTPARMS (locations 351 to 550) in 
the section of this manual describing the inputs. 

This program will generate the QWMU vector for sinusoidal gust 
described in Reference 1, pp 88-90. It requires these inputs: 

J, advance ratio 

NCP, number of chordwise panels (LSTPARMS, location 2) 

NSM, number of spanwise modes (LSTPARMS, location 4) 

QZAR, spanwise locations of control point radii (LSTPARMS, 
locations 

BD, blade chord/diameter at the QZAR points 

QQ, order of unsteady loading harmonic (LSTPARMS, location 21) 

WZERO, peak amplitude of gust velocity divided by rotor tip speed 

This program will print the elements of the QWMU vector, in a 
form suitable for use in the LSTPARMS section of the input. 

The second listing (Figure 52) is the block data, CRPBLK, in 
which the default values for the items described under the 
CRPPARMS input section are set. This listing has the default 
value of CRPTOL (CRPPARMS, location 5) set to 0.01. The default 
values of the other variables can also be set to different values 
here. 
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Table I 

Common block locations 


Common Block 
Name 

Size in bytes 

subprogram {type} 

airalz 

2364 

airbk3 {block data} 
air24 {subroutine} 
airdrg {subroutine} 
airflx {subroutine} 
airoff {subroutine} 
airp02 {subroutine} 
casarf {subroutine} 
drag24 {subroutine} 
isoafl {subroutine} 
isoarf {subroutine} 
lift24 {subroutine} 
swparf {subroutine} 
zeroal {subroutine} 

aircdf 

312 

airbk2 {block data} 
air24 {subroutine}} 
airdrg {subroutine} 
airflx {subroutine} 
airoff {subroutine} 
airp02 {subroutine} 
casarf {subroutine} 
drag24 {subroutine} 
isoafl {subroutine} 
isoarf {subroutine} 
lift24 {subroutine} 
swparf {subroutine} 
zeroal {subroutine} 

aircdm 

392 

airbk3 {block data} 
air24 {subroutine} 
airdrg {subroutine} 
airflx {subroutine} 
airoff {subroutine} 
airp02 {subroutine} 
casarf {subroutine} 
drag24 {subroutine} 
isoafl {subroutine} 
isoarf {subroutine} 
lif t24 {subroutine} 
swparf {subroutine} 
zeroal {subroutine} 
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Table I (continued) 


airdat 

4000 

air24 {subroutine} 
airdgr {subroutine} 
airflx {subroutine} 
airoff {subroutine} 
airp02 {subroutine} 
casarf {subroutine} 
drag24 {subroutine} 
isoafl {subroutine} 
isoarf {subroutine} 
lift24 {subroutine} 
swparf {subroutine} 
zeroal {subroutine} 

airpOl 

400 

airbkl {block data} 
airdrg {subroutine} 
airprm {subroutine} 
lstuns {subroutine} 

bdscOl 

644 

bds018 {subroutine} 
bldbk3 {block data} 

bdsc02 

8 

bds018 {subroutine} 

bdsc03 

3628 

bds014 {subroutine} 
bldbk4 {block data} 

besar 

6400000 

bessav {subroutine} 
bsint {subroutine} 
bsjint {subroutine} 
bsy int { subrout ine } 

besdel 

8 

bessav {subroutine} 
bsint {subroutine} 
bsjint {subroutine} 
bsyint {subroutine} 

bldgOl 

4000 

bldbk {block data} 
bldctl {subroutine} 
bldgil {subroutine} 
fetch {subroutine} 
retrev {function} 
store {subroutine} 

bldg02 

8000 

fetch {subroutine} 
retrev {function} 
store {subroutine} 



Table I (continued) 


bldg03 

60800 

bldct3 (subroutine) 
bldgi3 (subroutine) 
fetch (subroutine) 
retrev (function) 
store (subroutine) 

bldgc3 

121600 

fetch (subroutine) 
store (subroutine) 

bldgeo 

24000 

bldbk2 (block data) 
bldctl (subroutine) 
bldct3 (subroutine) 
bldgil (subroutine) 
bldgi2 (subroutine) 
bldgi3 (subroutine) 
bldgpl (subroutine) 
fetch (subroutine) 
lstuns (subroutine) 
mcafa (subroutine) 
retrev (function) 
store (subroutine) 

bldgxo 

48000 

fetch (subroutine) 
retrev (function) 
store (subroutine) 

com002 

72 

ggsOll (subroutine) 
ggsl07 (subroutine) 

crppOl 

36 

crpblk (block data) 
crpprm (subroutine) 
lstuns (subroutine) 
lstwc (subroutine) 
newpgl (subroutine) 
nrkOOl (subroutine) 
retrev (function) 
wakOOl (subroutine) 

crpsrl 

24 

velgrd (subroutine) 

crpsr2 

4 

velgrd (subroutine) 
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Table I (continued) 


dtetme 


32 


headrl 


728 


headr2 


4 


hsgspl 


4424 


intpOl 


64 


kclsav 


40 


airprm {subroutine} 
bldgem {subroutine} 
bldgil {subroutine} 
bldgi2 {subroutine} 
bldgi3 {subroutine} 
bldgpl {subroutine} 
hdrbkl {block data} 
header {subroutine} 
inptrw {subroutine} 
intfrl {subroutine} 
intprm {subroutine} 
lstctl {subroutine} 
lstprm {subroutine} 
lstuns {subroutine} 
MAIN {main routine} 
newpgl {subroutine} 
noisez {subroutine} 
runprm {subroutine} 
vtxprm {subroutine} 
wakOOl {subroutine} 
wakprl {subroutine} 
wakprm {subroutine] 

hdrbk2 {block data} 
header {subroutine} 
lstctl {subroutine} 
lstuns {subroutine} 
newpgl {subroutine} 

hdrbk2 {block data} 
header {subroutine} 
lstctl {subroutine} 
lstuns {subroutine} 
newpgl {subroutine} 

lstswl {subroutine} 
noznhf {subroutine} 
qsplin {subroutine} 
vtxinp {subroutine} 

intblk {block data} 
intfrs {subroutine} 
intfrl {subroutine} 
lstctl {subroutine} 
lstuns {subroutine} 

isoarf {subroutine} 


71 


Table I (continued) 


IstpOl 

2800 

fetch {subroutine) 
lstbkl {block data) 
lstprm {subroutine} 
lstuns {subroutine} 
retrev {function} 
store {subroutine} 

lstp02 

5600 

fetch {subroutine} 
retrev {function} 
store {subroutine} 

IstrOl 

29800 

airdrg {subroutine} 
fetch {subroutine} 
intrfl {subroutine} 
lstuns {subroutine} 
lstwc {subroutine} 
store {subroutine} 

lstr02 

4400 

fetch {subroutine} 
intfrl {subroutine} 
lstuns {subroutine} 
store {subroutine} 
wakOOl {subroutine} 

lstsr2 

8800 

fetch {subroutine} 
store {subroutine} 

lstspr 

59600 

fetch { subrout ine } 
store { subrout ine } 

nozdat 

632 

nozbkl {block data} 
nozclc {subroutine} 
nozprm {subroutine} 

nozsl6 

6408 

nozbk2 {block data} 
nozffc {subroutine} 
noznfc {subroutine} 

output 

4 

bds014 {subroutine} 

provrs 

4 

idblk {block data} 
lstuns {subroutine} 
newpgl {subroutine} 

rncpOl 

400 

lstctl {subroutine} 
lstuns {subroutine} 
mcafa {subroutine} 
runblk {block data} 
runprm {subroutine} 
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Table I (continued) 


vliicl 

4900 

lstctl {subroutine} 
Is tuns {subroutine} 
vel002 {subroutine} 
velbkl {block data} 
velgrd {subroutine} 
vttbin {subroutine} 

vtxc02 

480 

fetch {subroutine} 
retrev {function} 
store {subroutine} 

vtxcom 

240 

fetch {subroutine} 
retrev {function} 
store {subroutine} 
vtxblk {block data} 
vtxout {subroutine} 
vtxprm {subroutine} 

wakpOl 

4000 

fetch {subroutine} 
lstuns {subroutine} 
retrev {function} 
store {subroutine} 
wakOOl {subroutine} 
wakblk {block data} 
wakprm {subroutine} 

wakp02 

8000 

fetch {subroutine} 
retrev {function} 
store {subroutine} 

wakrOl 

24000 

fetch {subroutine} 
intfrl {subroutine} 
lstuns {subroutine} 
store {subroutine} 
wakOOl {subroutine} 
wakpr 1 { subrout ine } 

vaksp2 

48000 

fetch {subroutine} 
store {subroutine} 

work 

2119852 

f271 {subroutine} 
f271no {subroutine} 
intfr3 {subroutine} 
lstctl {subroutine} 
lstuns {subroutine} 
nrkOOl {subroutine} 
wak003 {subroutine} 

workl 

3200 

f 27 1 {subroutine} 
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Table Z (continued) 


worksp 

xindex 


100000 

12 


f271no {subroutine} 
lstuns {subroutine} 

MAIN {main routine} 

bds014 {subroutine] 
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Table II 

Locations of IMSL routines used within UCAP 


IHSL Routine 

Used by 



bsio 

ikO, iksubO, stdiki, sum3. 

suin4 , tikO 

bsil 

sum3 



bsins 

nozhnk, wake 



bsjo 

iksubO 



bsjl 

wingf 



bsjns 

bessav, bsint, 

bsjint, nozffc, nozhnk 

bskO 

ikO, iksubO, stdiki, sum3. 

sum4, tikO 

bskl 

sum3 , wingf 



bskes 

nozhnk, wake 



bsys 

bessav, bsint. 

bsyint, iksubO, nozhnk 

csint 

fOOl 



csval 

fOOl, ffun 



ctime 

jobtim 



dtiroe 

MAIN, endjob 



erfc 

imnsub, imnsin 



fftcb 

fftcx, fftcxO, 

f ftcxn 


iidex 

fetch, retrev, 

store 


iicsr 

fetch, retrev. 

store 


iwkin 

MAIN 



lincg 

kmatrx 



linrg 

f271no, kmatrO 



tdate 

MAIN 




noznfc 


wingf 
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Table III 

Locations where AIMAG is used within UCAP 


ansnl 

bds024 

bound 

f ftcxn 

iksub 

iksubO 

lstctl 

nozf fc 

noznj2 

noznfc 

noznj2 

no z out 

nozstt 

nozunt 

r4aray 

wingf 
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V xp oti AV X /V 0 due to potential loading of forward rotor. 

Vtpoti Vj/Vo due to potential loading of forward rotor. 

V ts1 V t /Vo due to secondary (non-potential) loading of forward rotor. 

V xs1 AV x /Vo due to secondary (non-potential) loading of forward rotor. 

Vx 2 i AV x A/ 0 from rear rotor loading at forward rotor. 

V b i AV x /Vo due to presence of center body. 

V 0 Freestream velocity. Normalized freestream velocity is unity. 

ayZ (jt/Ji)»(r/Ri), where J is advance ratio, is tip radius, r is local radius. 

<t> Advance angle, measured from plane of rotation. 

© t1 Camber plus twist angle for each panel, measured from plane of rotation. 

Wt Turning due to blade thickness. 

Wri Turning required for boundary condition (no flow through the blade) to be met. 
AV X Change in axial velocity 

V t Tangential (swirl) velocity 


Figure 1 (part a) Forward Rotor Velocity Diagram 


DG1217001ag 
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v xpot 2 AV X /V 0 ,j ue t 0 potential loading of rear rotor. 
v tpot 2 V, /V 0 due to potential loading of rear rotor. 

v ts 2 V,/V 0 due to secondary (non-potential) loading of rear rotor. 
v xs2 AV X /V 0 due to secondary (non-potential) loading of rear rotor. 
v xi2 AV X /V 0 from forward rotor loading at rear rotor. 

V,i 2 V x /V 0 from forward rotor loading at rear rotor. 

V b2 AV X /V 0 due to presence of center body. 

V 0 Freestream velocity. Normalized freestream velocity is unity. 

a 2 -Z (7r/J 2 ) , (r/R2), where J is advance ratio, R 2 is tip radius, r is local radius. 

Advance angle, measured from plane of rotation. 

© t 2 Camber plus twist angle for each panel, measured from plane of rotation. 

W t Turning due to blade thickness. 

Wp Turning required for boundary condition (no flow through the blade) to be met. 
AV X Change in axial velocity 

Vt Tangential (swirl) velocity 

Figure 1 (part b) Rear Rotor Velocity Diagram 


I 
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A 

B 

C 

D 

E 

F 

G 

H 

J 

K 

L 

Steps 

M 

N 

P 

Q 

R 

S 

T 

U 

V 

W 

X 


Calculate axisymmetric streamlines. This is done off line. 

Calculate, off line, streamline angles where streamlines cross blade 
mid-chord line. These angles are measured from the axis of rotation. 

Calculate table of camber plus twist angles, along conical surfaces 
defined by streamline angles calculated in step B, for front and rear 
rotors . 

Calculate the location of the rear rotor control points in relation to 
stream tubes defined by the axisymmetric streamlines. 

Calculate forward (or single) rotor influence coefficient matrix. Store 
for later use. 

Calculate turning angle due to blade thickness for forward (or single) 
rotor. Store for later use. 

Calculate potential loading. L=[K]'xW, where [K)' 1 is influence coeffi- 
cient matrix, and W is the turning angle required to satisfy the 
boundary conditions. 

Calculate induced velocity from potential loading using momentum theory. 
Calculate total loading, including vortex effects. 

Calculate induced velocity from total loading, including vortex effects. 

Calculate axial AV/V® and swirl, Vfo/ v o. at rear rotor mapping points, 
calculated in step D, above. These velocities will be applied at the 
rear rotor control points, as shown in Figure 1. 


if through V are only applied In counter-rotation . For single rotation, 

jump to step W. 

Calculate locations of forward rotor control points relative to annular 
stream tubes defined by the axisymmetric streamlines. 

Calculate rear rotor influence coefficient matrix. Store it for later 
use. 

Calculate turning angle due to blade thickness. 

Calculate potential loading, as in step G, above. 

Calculate induced velocity from potential loading, using momentum theory. 
Calculate total loading, including vortex effects. 

Calculate induced velocity field from total loading, including vortex 
effects . 

Calculate axial AV/V. induced by rear rotor at forward rotor mapping 
points, calculated in step M, above. Note there is no swirl ahead of a 
rotor. 

Calculate synthetic streamline radial locations in order to enforce 
incompressible continuity, V*V=0. 

Convergence monitor: converged if | c pi -c pH | is less than a limiting value 

for both rotors. 

Noise module: acoustic analysis 

Figure 2: Conceptual flow chart of UCAP 
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Figure 3A: UCAP Calling Tree 
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Figure 3A: UCAP Calling Tree (continued) 
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Figure 3A: UCAP Calling Tree (continued) 
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Figure 3A: UCAP Calling Tree (continued) 
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Figure 3A: UCAP Calling Tree (continued) 
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Figure 3A: UCAP Calling Tree (continued) 
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Figure 3A: UCAP Calling Tree (continued) 
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9 


Oft AG 2 4 


A 


LNEARI 


5 


CASARF 


4 


ISOARF 


7 


ISOAFL 


8 


AIR24 


9 


LIFT24 


A 


LNEARI 


9 


DRAG24 


A 


LNEARI 


8 


AIROFF 


9 


AIR24 


A 


LIFT24 

/ 

8 


LNEARI 


A 


DR AO 24 


8 


LNEARI 
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4 

STARC 



7 

UN I NT 



4 

zeroal 



7 

I SOAR F 



0 

ISOAFl 



0 

AIR24 



A 

nrr24 



0 

LNEARl 



A 

DRAG24 



0 

LNEARl 



1 

AIROFF 



A 

AIR24 



0 

l IFT24 



C 

LNEARl 



0 

DRAG 24 



C 

LNEARl 



t 

DZROAL 



7 

UNBAR 



7 

UN I NT 



4 

CLFACT 



7 

UNBAR 



5 

LNEARl 



5 

NENPOl 



3 

AIRP02 



A 

PftNTIN 



4 

PRNTXX 



3 

PRNTXX 



l BLDGEH 




2 BLDG 1 1 



3 

INPTC1 



A 

END JOB 



5 

□TIME 



A 

INPTC2 



3 

NEWPOl 



3 

PRNTXX 



3 

BLDCTl 



A 

QSPLIN 



5 

SPFUNC 



5 

SPCOEF 



4 

SPL5D3 



5 

SPSVAL 



4 

SPL5EV 



5 

5PMVAL 



4 

SPL5EV 



A 

UVINT 



A 

NENPOl 



A 

PR NT IN 



4 

PRNTXX 



A 

BDS052 



5 

EXIT 



5 

BDS004 



4 

BDS01S 



7 

BDSQ 1 4 



0 

BDSOIO 



0 

CAMTHA 



A 

EXIT 



0 

BDSO10 



0 

HS1 



0 

EXIT 



0 

PF1 



0 

EXIT 



0 

MS 2 



0 

EXIT 



0 

H455S 



0 

EXIT 



0 

PF2 



0 

EXIT 



6 

BDSOll 



0 

BDS020 



0 

BDSQ 1 0 



0 

BDSQ 13 



0 

MLINE 



0 

BDSQ 1 0 



0 

BDSQ 13 



0 

CAMTHA 



0 

EXIT 



0 

BDS012 



0 

BDSQ 24 



0 

BDSOlft 



A 

ERRSXT 



0 

BDS025 



A 

BDS024 



A 

BDSQ23 



0 

BDSQ 20 



7 

BDS024 



0 

BDSOlft 



0 

ERRSXT 



0 

BDS025 



0 

BDSQ 2ft 



0 

BDSQ 23 



4 

BDS104 



7 

BDS0S4 



0 

BDSO04 



0 

BDSO05 



0 

BDSO07 



0 

BDSO00 



0 

BDSDB0 



0 

BDSO0S 



5 

CCS 107 



5 

GGSOll 



5 

BDS024 



4 

BDSOlft 



7 

ERRSXT 



4 

BDS025 



7 

BDS024 



7 

BDSQ 23 



5 

BDS053 



4 

EDS 024 



7 

BDSOlft 



0 

ERRSXT 



7 

BDS025 



0 

BDS024 



B 

BDS023 



5 

BDSOlft 



4 

ERRSXT 



5 

GGS252 



4 

0GS253 



4 

GGS25A 



5 

BDS05A 



4 

EXIT 



A 

UNINT 



A 

BLDOCA 
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s 

3 

4 

4 


5 SLOCUB 

5 NEWPOl 

BIDOPI 

4 NEWPG1 

4 PRNTXX 

KD0I2 

INPTCl 

4 ENDJOB 

S DTIME 

4 INPTC2 

NEWPOl 
BLDGP1 

NEWPG1 
PRNTXX 
UN I NT 
BUXJI5 

3 INPTCl 

4 ENDJOB 

5 DTIME 

4 INPTC2 

3 NEWPGl 

3 PRNTXX 

3 EXIT 

3 BLDCT3 

4 PRNTIN 

4 BDSA52 

5 BDS024 

4 BDSO It 

7 ERRSXT 

4 BDS025 

7 BDS024 

7 BDS023 

5 BDS053 

4 BDS024 

7 BDSOlt 

a ERRSXT 

7 BDS02S 

a BDS024 

a BDSO 23 

5 BDSOlt 

4 ERRSXT 

5 GGS252 

4 GGS253 

4 GGS254 

5 BDS054 

EXIT 
MEWPG1 
PRNTXX 
UNINT 
BLDOCA 

SIOCUB 
NEHPG1 

1LDOP1 

NEWPG1 
PRNTXX 

STORE 
CRPPRH 

INPTCl 

ENDJOB 
DTIME 
rNPTC2 


4 

4 

4 

4 

4 

5 

5 

3 

4 

4 


2 NEWPOl 

2 PRNTXX 

1 ENDCAS 

2 DTIME 

1 ENDJOB 

2 DTIME 

1 HEADER 

2 NEWPGl 

1 INTPRM 

2 INPTCl 

3 ENDJOB 

4 DTIME 

3 INPTC2 

2 NEWPOl 

2 PRNTXX 

1 LSTPRM 

2 FETCH 

2 INPTCl 

3 ENDJOB 

4 DTIME 

3 INPTC2 

2 NEWPGl 

2 PRNTXX 

2 PRNTCS 

2 STORE 

I NO I SEX 

2 WOZCLC 

3 NEWPOl 

3 JOBTIM 

3 MOZINP 

3 NOZCHK 

4 NOZNJ2 

5 NOZTCM 

4 WOZXKC 

5 MNX 

5 NO 2 AMP 

4 W02HNK 

7 DBS I NS 

7 DBSKES 

7 DBSJNS 

7 DBSYS 

3 NOZFFC 

4 NOZTRP 

5 EXIT 

4 DBSJNS 

4 NOZTCM 

4 NOZAIR 

4 NOZMPL. 

5 NOZM4F 

4 OSPLIN 

7 SPFUNC 

7 SPCOEF 

a S PL SB 3 

7 SPSVAL 

a spisev 

7 SPMVAl 

a SPL5EV 

4 NOZHFL 

4 N02TRL 

4 NOZBLC 

4 NOZOPL 


Figure 3A: UCAP Calling Tree (continued) 
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6 

MOZTRL 


4 

NOZUNT 


4 

MOZSTT 


3 

NOZNFC 


4 

MOZTRP 


5 

EXIT 

_ _ 

4 

MOZHNK 

s 

5 

DBS IMS 


5 

DBSXES 


5 

DBSJMS 


5 

DBSVS 


4 

MOTTO 


4 

NOZAIR 


4 

NOZHPL 


5 

MOZNHF 


4 

QSPLIM 


7 

SPFUNC 


7 

SPCOEF 


a 

SPL5D3 


7 

SPSVAL 


a 

SPL5EV 


7 

SPMVAl 


6 

SPL5EV 


4 

NOZHFl 


4 

NOZTRL 


4 

MOZBLC 


4 

MOZQPL 

' 

4 

DBSJMS 


4 

MOZTRL 


4 

MOZUMT 



MOZSTT 


3 

NOZOUT 


4 

MEWPGI 


4 

MOZOU2 


5 

MEWPGI 


5 

MOZUMT 

P. ^ 

5 

MOZSTT 


5 

NOZOU3 


i * 

MEWPGI 


4 

MOZUMT 


1 NOZPRM 

-— 

2 

IMPTCI 


3 

EMDJOB 

- - 

4 

DTIME 


3 

INPTC2 


2 

MEW POL 


2 

PRMTXX 


I RUNPRM 


2 

IMPTCI 


3 

EMDJOB 


4 

DTIME 

"" "" 

3 

IMPTC2 


2 

MEWPGI 


2 

PRMTXX 


l VELGRD 


2 

IMPTCI 

_ _ 

3 

EMDJOB 

’ : 

4 

DTIME 

— — 

3 

IMPTC2 


2 

MEWPGI 

sir 

2 

VEL002 


3 

PRMTXX 


1 VTXPRH 


2 

FETCH 

■ — 

2 

MEWPGI 


2 

INPTC1 


3 

ENDJOB 


4 

DTIME 


3 

IMPTC2 

— 

2 

PRMTXX 


2 

STORE 

1...Z 

1 WAX001 


2 

FETCH 


2 

WAX 002 


3 

ERRSVT 


5 

MEWPGI 


3 

WAKPR1 

m-. -a 

4 

FETCH 

“ - 

4 

PRWTDP 


4 

PRMTIM 


3 

EXIT 


3 

F001 


4 

DWORM 


5 

DMMX 



4 

IFXPLT 


4 

DCS I WT 


4 

IMILU 

Z=T-~1 

3 

WAK003 


4 

EXIT 


4 

MM SUB 


5 

PI 


5 

P2 

_ _ 

5 

CHECXM 


4 

NRK001 


5 

JOBTIM 


5 

MEWPGI 


5 

DTABLP 


5 

PRNTCD 


5 

BESSAV 



4 

EXIT 

- 

4 

DBSJMS 

t— ^ 

4 

DBSVS 

tjj2§ 

5 

BLEND 


4 

DUNLIN 


4 

IMILU 


5 

MUNUAR 


5 

WRK002 

-— = 

4 

NRK007 


7 

EXIT 


1 

NRX004 


4 

MAX 0 05 


4 

MRK003 


7 

ZODEF 


7 

DUMLIM 

-- 3 

7 

MRK009 

* * 

6 

1X0 


6 

FFTCXN 


1 

DFFTCB 


B 

PLTFFO 


f 

IMILU 


4 

JOBTIM 


4 

MIK004 

r=1 L 
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ZOOEF 

DUWLIM 

MRKOOa 

IKMCW 

RSJIWT 

OBSJMS 

BSVINT 

DBSYS 

JVNEW 

IKMHP 

JVCOMB 

FFTCHN 

DFFTCB 

PLTFFO 

IHILU 

CHECK 



Figure 3A: UCAP Calling Tree (continued) 
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ROUND 0 IS CALLED BY 
BSINT IS CALLED BY 
BSJINT IS CALLED BY 
BSYINT IS CALLED BY 
CAHTKA IS CALLED BY 
CASARF IS CALLED BY 
CENTER IS CALLED BY 
CHECX IS CALLED BY 
CHECKM IS CALLED BY 
CKOUT IS CALLED BY 
CLFACT ISj CALLED BY 
CONVRO IT CALLED BY 
CPMAP IS CALLED BY 
CPYSDP IS CALLED BY 
CROSSP IS CALLED BY 
CRPPRH IS CALLED BY 
CSPSUB IS CALLED BY 
CTRAP IS CALLED BY 
CUNLIN IS CALLED BY 
DEF20 IS CALLED BY 
DMNX IS CALLED BY 
DNORM IS CALLED BY 
DRAG24 IS CALLED BY 
DSLOCB IS CALLED BY 
0TABLP IS CALLED BY 
DTRAP IS CALLED BY 
DUNBAR IS CALLED BY 
OUNLIN IS CALLED BY 
DUNLIN IS CALLED BY 
DUNLIN IS CALLED BY 
DZROAL IS CALLED BY 
ENDCAS IS CALLED BY 
ENDJOB IS CALLED BY 
ENDPT IS CALLED BY 
ERRSXT IS CALLED BY 
ERRSVT IS CALLED BY 
ETACAL IS CALLED BY 
FETCH IS CALLED BY 
FEXP 15 CALLED BY 
FFTCX IS CALLED BY 
FFTCXN IS CALLED BY 
FFTCX 0 IS CALLED BY 
FFUN IS CALLED BY 
F001 IS CALLED BY 
F27I IS CALLED BY 
F271NO IS CALLED BY 
QETXMB IS CALLED BY 
GGS011 IS CALLED BY 
005107 IS CALLED BY 
GGS252 IS CALLED BY 
GGS255 IS CALLED BY 
COS 25 A IS CALLED BY 
HEADER IS CALLED BY 
HNSUB IS CALLED BY 
H$1 IS CALLED BY 
HS2 IS CALLED BY 
H455S IS CALLED BY 
IBIN IS CALLED BY 
IFXPLT IS CALLED BY 
IIEXPS IS CALLED BY 
IIEXW IS CALLED BY : 


: KHATRO 
: CSPSUB 
i IKNEW 
i I KNEW 

l BDSOIO BDSOU 
AIRFLX SWPARF 
BOUND 

HSUH NRKOOI THICKS WAKE 

HHSUB 

INTLIN 

CASARF 


WAKEOO WAKESM HAKES 0 


ANSO AWSOMO 

ANSN1 ANSONO ANSON 1 BLEND BOUND BO UNDO CSPSUB F271 INTERS 
KHATRO LSTWC NRK003 NRKOOi THICK 0 THICKS HAKE WAKEQ9 HAKESN 
HAKHO HAKHOl HAK004 


LSTDY LSTPRM LSTUNS HCAFA VTXPRH HAKPRM HAKPRI WAK001 

IWFFNS HINOIN 
NRKMB NRKOOI 
INF INFFNO 

INTFR2 HAX002 

LSTDY LSTUNS 

LSTDY 

LSTCTL 

BDS052 

BDS052 

BDSA52 BDS052 
QOS 252 
0CS252 
HAIM 

INTFRS HAXOOS 
BDS014 
BDSOU 
BDSOU 

F001 NRK01A NRKOIO 
BO UNDO THICKO THICKS 
BOUND 


Figure 3B: UCAP Subroutine Cross Reference Listing 
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IS CALLED 

' IS CALLED 

I IS called 
IS CALLED 

, IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS called 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 

! IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 
IS CALLED 


BOUND 

QSPLIN 

IKSUSO 

QSPLIN 

SPCOEF 

SPMVAL SPSVAL 

QSPLIN 

QSPLIN 

LSTCTL 

CAS ARE 


rldgeh lstdv lstprh LSTUNS MCAFA vtxprm WAXPRH WAX 0 0 1 


DZROAL 

BLDCT3 BLDGI2 DZROAL STARC 

ANSO ANSONO BLDCT1 SLFWD VTTBIN 


MAIN 

VELGRO 

LSTWC 

LSTCTL 

VTXCAL 

LSTDY 

VTXCAL 

VTXCAL 

VTXCAL 

VTXCAL 

VTXCAL 

VTXCAL 

LSTDV LSTUNS 

LSTDV VTXCAL 

VTXCAL 

MAIN 

LSTDV 

KHATRX 

KMATRO 

WAKE 

WAXEQO 

KMATRO 

WAXMO 

MAIN 

WAXQ02 

MAIN 

WAX 0 0 1 

WAX 0 02 

INTFR3 WAXOOS 
BOUND 


WINOF IS CALLED BY : WING IN 
WINGIN IS CALLED BY : WINO 
2EROAL IS CALLED BV 1 CASARF 

20DEF IS CALLED BY : BOUNDS NRX803 NIRKO0* THICXO THICKS 


Figure 3B: UCAP Subroutine Cross Reference Listing (continued) 
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■ *** TSO FOREGROUND HARDCOPY ■■■» 10:34:55 33112 



DSNAME-TSOC04 I . XREF2 . DATA 




SUBROUTINE REFERENCE ANO 

PURPOSE LISTING 



SUBROUTINE ABSVAl 
SUBROUTINE ADJUST 

PURPOSE : 

PURPOSE : 

CALL <S) : NEWPOl UVINT 



SUBROUTINE AIRDRO 

PURPOSE l 

CALL (SI : AIRFLX NEWPOl AIRP02 PRNTXX 



SUBROUTINE AIRFLX 

PURPOSE : 

CALL (S3 : ISOAftF CASARF SWPARF 



SUBROUTINE AIROFF 
SUBROUTINE AIRPRM 

PURPOSE l 
CALL (S3 : AIR24 
PURPOSE i 

CALL (S3 : INPTC1 NEWPOl PRNTXX AIRDRO 



SUBROUTINE AIRP02 

PURPOSE : 

CALL (S) : PRNTIN PRNTXX 



SUBROUTINE AIR24 

PURPOSE l 

CALL (S) : LIFT24 DRAG24 



SUBROUTINE ANS 

PURPOSE t 

CALL (SI : LSTWC DUNBAR DSLOCB DUNLIN ANSN1 



SUBROUTINE ANSN1 

PURPOSE t 

CALL (S3 : NEWPOl DUNLIN INILU 



SUBROUTINE ANSO 

PURPOSE : 

CALL (S3 : LSTWC DUNBAR ANS0N1 NEWPOl UVINT 

AIRDRO 


SUBROUTINE ANSONO 

PURPOSE : 

CALL (S3 : DUNLIN DUNBAR AWS0N1 NCWPGl UVINT 

AIRDRO 


SUBROUTINE ANS 9N1 

SUBROUTINE ARCTAN 
SUBROUTINE A1K2 
SUBROUTINE BDSA52 

PURPOSE l 

CALL (S3 : DUNLIN IMILU 
PURPOSE » 

PURPOSE : 

PURPOSE 1 



CALL (S) : BDS024 BDS053 BDSOli GGS252 BDSC54 


SUBROUTINE BDS004 

PURPOSE J 

CALL ($> : BDS015 BDS104 



SUBROUTINE BDS010 

PURPOSE i 

CALL (S3 : CAMTHA 



SUBROUTINE BDSOli 

PURPOSE t 
CALL (S3 ; BDS820 



SUBROUTINE BDS012 

PURPOSE 1 



SUBROUTINE BOSS 13 

PURPOSE : 



SUBROUTINE BDS014 

PURPOSE ; 

CALL (S3 : HLINE CAMTHA BDS012 BDS824 BDS020 
CALL (S3 : BDS010 BDS01B HS1 PF1 HS2 

H455S PF2 

BDSOM BDSOli BDS0I3 

SUBROUTINE BDS015 

PURPOSE : 

CALL (S3 : BDSQ14 BDS024 



SUBROUTINE BDSOli 

PURPOSE : 

CALL (S3 : ERRSXT 



SUBROUTINE 8DS01B 

PURPOSE : 



SUBROUTINE BDSOli 

PURPOSE : 



SUBROUTINE BDS020 

PURPOSE : 



SUBROUTINE BDS02S 

PURPOSE : 



SUBROUTINE BDS024 

PURPOSE : 

CALL (S3 : BDS014 BDS025 



SUBROUTINE BDS025 

PURPOSE : 

CALL (S3 : BDS024 BDS023 



SUBROUTINE BDS924 

PURPOSE : 



SUBROUTINE BDS052 

PURPOSE : 




CALL (S3 : EXIT BDS004 GOS107 OOSOll BDS024 

BDS053 BDS014 

GGS252 BDS054 

SUBROUTINE BDS053 

PURPOSE : 

CALI (S3 ; 80S024 



SUBROUTINE BDS054 

PURPOSE t 



SUBROUTINE BDS0A6 

PURPOSE : 

CALL (S3 : BDSOli BDS015 BDS0B1 BDSOBB 



SUBROUTINE BDSOBB 

PURPOSE : 



SUBROUTINE BDSOdf 

PURPOSE : 

CALL (S3 : BDSOBB 



SUBROUTINE BDS015 

PURPOSE i 
CALL (S3 : BDS007 



SUBROUTINE BDS014 

PURPOSE : 



SUBROUTINE BDS017 

PURPOSE : 



SUBROUTINE BDS104 

PURPOSE : 

CALL (S3 : BDSQ64 



SUBROUTINE BESSAV 

PURPOSE : 

CALL (S3 : EXIT DBSJNS DBS VS 



SUBROUTINE BLDCT1 

PURPOSE : 

CALL (S3 : QSPLIN UVINT NEWPOl PRNTIN PRNTXX 

BDS052 EXIT 

UN I NT BLDOCA 

SUBROUTINE BLDCT3 

PURPOSE : 

CALL (S3 : PRNTIN BDSA52 EXIT NEWPOl PRNTXX 

UN I NT BLDOCA 


SUBROUTINE BLDOCA 

PURPOSE i 

CALL (S3 : SLOCUB NEWPOl 



SUBROUTINE BLDOEM 

PURPOSE J 

CALL (S3 : BLDOI1 BLDOI2 BLDOI3 STORE 



SUBROUTINE BLDOI1 

PURPOSE : 

CALL (S3 : INPTC1 NEWPOl PRNTXX BLDCT1 BLDCPI 



SUBROUTINE 8LDOI2 

PURPOSE : 

CALL (S3 : INPTCl NEWPOl BLDOPI UNIMT 



SUBROUTINE BLDG I 3 

PURPOSE f 

CALL (S3 : INPTCl NEWPOl PRNTXX EXIT BLDCT3 

BLDGPl 


SUBROUTINE BL DO PI 

PURPOSE : 

CALL (S3 : NEWPOl PRNTXX 



SUBROUTINE BLEND 

PURPOSE : 

CALL (S3 : DUNLIN IMILU 



SUBROUTINE BOUND 

PURPOSE J 

CALL (S3 : SMWAJCO SHWAK DUNLIN IIEXYY DCFZO 

HSUH Wt WO 

CENTER RAD INT 

SUBROUTINE BOUNDO 

PURPOSE : 

CALL (S3 : DUNLIN IIEXPS ZODEF INFFNO 



SUBROUTINE BSINT 

PURPOSE t 

CALL (S3 : DBSJNS DBSYS 



SUBROUTINE BSJINT 

PURPOSE i 

CALL (S3 : DBSJNS 



SUBROUTINE B5YIMT 

PURPOSE l 
CALL (S3 : DBSYS 



SUBROUTINE CAMTHA 

PURPOSE : 

CALL (S) i EXIT 



SUBROUTINE CASaRF 

PURPOSE : 

CALL (S3 s ISOARF STARC ZEROAL DZROAL CLFACT 



SUBROUTINE CENTER 

PURPOSE S 



SUBROUTINE CHECX 

PURPOSE : 



SUBROUTINE CHECKH 

PURPOSE : 


- 

SUBROUTINE CKOUT 

PURPOSE : 



SUBROUTINE CLFACT 

PURPOSE : 

CALL (S) : UNBAR 



SUBROUTINE CONVRG 

PURPOSE : 



SUBROUTINE CPHAP 

PURPOSE i 



SUBROUTINE CPYSDP 

PURPOSE : 



SUBROUTINE crossp 

PURPOSE : 



SUBROUTINE CRPPRH 

PURPOSE i 

CALL (S) i INPTCl NEWPOl PRNTXX 






Figure 3C: 
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SUBROUTINE LSTDY 


SUBROUTINE LSTPRM 

SUBROUTINE LSTSW1 

SUBROUTINE L STUNS 

SUBROUTINE LSTWC 

SUBROUTINE LSTXY2 
SUBROUTINE MAIN 


SUBROUTINE HCAFA 

SUBROUTINE MEAN 
SUBROUTINE HLINE 

SUBROUTINE MNX 
SUBROUTINE MOMCAl 
SUBROUTINE MSUH 

SUBROUTINE MUNUAR 
SUBROUTINE NEWPOl 
SUBROUTINE NO I SEX 

SUBROUTINE NOZAIR 
SUBROUTINE WOZAHP 

SUBROUTINE NOZBIC 
SUBROUTINE N02CHK 

SUBROUTINE N02CLC 

SUBROUTINE MOZFFC 

SUBROUTINE NOZHFl 
SUBROUTINE N02HNK 

SUBROUTINE N02INP 
SUBROUTINE NOZKKC 

SUBROUTINE NOZHPL 

SUBROUTINE HOZNFC 


SUBROUTINE NOZNHF 
SUBROUTINE NOZNJ2 
SUBROUTINE WOZOUT 
SUBROUTINE NOZOU2 
SUBROUTINE M0Z0U3 


PURPOSE i 

CAU (S) J INPTCI NEWPOl PRNTXX 
PURPOSE i 
PURPOSE t 
PURPOSE : 

PURPOSE : 

PURPOSE i 
CALL (S) : EXIT 
PURPOSE I 
PURPOSE ; 

CALL (S) : JOBTIH NEWPOl DTABLP PRNTCD BESSAV BLEND MUNUAR NRK002 
PURPOSE : 

CALL (S) t NRK007 NRK0O4 NRK005 NRK003 JOBTIH NRKOQt 
PURPOSE ( 

CALL CS) i ZODCF DUNLIN NRKQIf 
PURPOSE : 

PURPOSE i 
PURPOSE t 

CALL CS) : 20DEF DUNLIN NRKOOB CHECX 
PURPOSE J 
CALL (S) i EXIT 
PURPOSE 5 

CALL (5) J I KNEW IKNWP FFTCXN PLTFFO 
PURPOSE i 

CALL CS) ; IK 0 FFTCXN PLTFFO 
PURPOSE J 

CALL <$) * NRK01B NEWPOl IHILU IFXPLT 
PURPOSE r 
PURPOSE t 

CALL CS) : EXIT NRKOl I NEWPOl IMILU IFXPLT 
PURPOSE : 

PURPOSE i 
CALL CS) : EXIT 
PURPOSE : 

CALL CS) i EXIT 
PURPOSE : 

CALL (S) : IHILU 
PURPOSE i 
CALL CS) : IMILU 
PURPOSE i 
PURPOSE 
PURPOSE : 

PURPOSE i 
PURPOSE i 
PURPOSE : 

CALL IS) : NEWPOl 
PURPOSE 1 
PURPOSE J 
PURPOSE : 

PURPOSE i 

CALL IS) : SPFUNC SPCOEF 5P5VAL 5PHVAL 
PURPOSE : 

CALL (S) : EXIT CUNLJN CTRAP 
PURPOSE s 
PURPOSE t 
PURPOSE : 

PURP O SE : 

CALL IS) : INPTCI NEWPOl PRNTXX 
PURPOSE I 
PURPOSE J 
PURPOSE : 


SUBROUTINE WOZPRH 

SUBROUTINE N02QPL 
SUBROUTINE NOZSTT 
SUBROUTINE NOZTCM 
SUBROUTINE N02TRL 
SUBROUTINE NOZTRP 

SUBROUTINE M02UNT 
SUBROUTINE NRK001 

SUBROUTINE NRK0O2 

SUBROUTINE WRK003 

SUBROUTINE NRKO04 
SUBROUTINE NRKOOB 
SUBROUTINE NRKOOt 

SUBROUTINE NRK007 

SUBROUTINE NRKOOB 

SUBROUTINE NRKOOt 

SUBROUTINE NRK01A 

SUBROUTINE NRK01B 
SUBROUTINE NRK010 

SUBROUTINE MRK 011 
SUBROUTINE PF1 

SUBROUTINE PF2 

SUBROUTINE PLTFFT 

SUBROUTINE PLTFFO 

SUBROUTINE PRNTCD 
SUBROUTINE PRNTCS 
SUBROUTINE PRNTDP 
SUBROUTINE PR NT IN 
SUBROUTINE PRNTXX 
SUBROUTINE PRTPRF 

SUBROUTINE PSIK2 
SUBROUTINE PI 
SUBROUTINE P2 
SUBROUTINE OSPLIN 

SUBROUTINE RAD I NT 

SUBROUTINE RETREV 
SUBROUTINE RJF 
SUBROUTINE RLCODE 
SUBROUTINE RUNPRM 

SUBROUTINE RAARAY 
SUBROUTINE SET 10 
SUBROUTINE SETLMT 


CALL (S) : SP2ER0 SET 10 OETXMB KINFIL THVFIL LSTDY NEWPOl LSTUNS SLFWD SIAFT 

PURPOSE i 

CALL (S) : STORE 

CALL (S) J FETCH NEWPOl PRNTXX F271NO F271 VTXCAL VTXOUT VTXNWT VTXWRD MOHCAL 
PURPOSE i 

CALL CS) : FETCH INPTCI NEWPOl PRNTXX PRNTCS STORE 
PURPOSE t 

CALL CS) J OSPLIN LSTCOA LSTXY2 
PURPOSE * 

CALL (S) : FETCH NEWPOl PRNTXX F271 VTXNWT STORE 
PURPOSE ; 

CALL CS) J NEWPGl PRNTDP PRNTIN DUNLIN VEL003 PRNTXX 
PURPOSE l 
PURPOSE J 

CALL <S) : ENOCAS ENDJOB HEADER INTPRM LSTPRM NO I SEX NOZPRM RUNPRM VELGRD VTXPRM 
CALL CS) i WAK001 WAKPRM 

CALL (S) s IWKIN JOBTIH DTIME TDATE INPTRW STORE LSTCTL AIRPRM BLDOEM CRPPRM 
PURPOSE J 

CALL CS) : FETCH NEWPOl TOLSTS PRNTXX PRNTIN ISTSW1 STORE 
PURPOSE : 

PURPOSE : 

CALL CS) : BDSO 19 BDS013 
PURPOSE l 
PURPOSE * 

PURPOSE * 

CALL CS) : INFFNS CHECX 

PURPOSE j 

PURPOSE J 

PURPOSE J 

CALL CS) : M02CLC 

PURPOSE i 

PURPOSE : 

CALL CS) : NOZHNK 
PURPOSE I 
PURPOSE : 

CALL CS) : N0ZNJ2 NOZKKC 
PURPOSE « 

CALL CS) : NEWPGl JOBTIH NOZINP NOZCHK NOZFFC NOZNFC NOZOUT 
PURPOSE J 

CALL CS) : NOZTRP DBSJNS NOZTCM NOZAIR NOZHPL MOZBL.C NOZQPL MOZTRL NOZUNT NOZSTT 
PURPOSE l 
PURPOSE J 

CALL CS) : DBSINS DBSKES DBSJNS DBSYS 
PURPOSE : 

PURPOSE : 

CALL CS) : »4X NOZAHP 
PURPOSE : 

CALL CS) : NOZNHF 
PURPOSE J 
CALL CS) : NOZSTT 

CALL (S) : NOZTRP NOZHNK NOZTCM NOZAIR NOZHPL NOZBLC NOZQPL DBSJNS N02TRL NOZUNT 
PURPOSE i 

CALL CS) : QSPLIN NOZHFL NOZTRL 

PURPOSE i 

CALL (S) : NOZTCM 

PURPOSE i 

CALL CS) : NEWPOl N0Z0U2 
PURPOSE r 

CALL CS) : NEWPGl NOZUNT NOZSTT NOZOU3 
PURPOSE t 

CALL CS) : NEWPOl NOZUNT 


B 


Figure 3C: UCAP Subroutine Reference and Purpose Listing (cont) 
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■ TSO FOREGROUND HARDCOPY **■■* 10:3«:SS 83112 
OSNAME-TSOOO*! .XREF3. DATA 

LABELED COMMON AREA CROSS REFERENCE LISTING 
[abel SUBROUTINE 



IS 

USED 

BY 

MAIN 





AIRALZ 

IS 

USED 

BY 

AIRDRG 

AIRFLX 

AIROFF 

AIRP02 

AIR2* 

AIRALZ 

IS 

USED 

BY 

SWPARF 

ZEROAL 




AIRCDF 

IS 

USED 

BY 

AIRDRG 

AIRFLX 

AIROFF 

AIRP02 

AIR2* 

AIRCDF 

IS 

USED 

BY 

SWPARF 

ZEROAL 




AIRCDM 

IS 

USED 

BY 

SWPARF 

ZERGAL 




AIRCDM 

IS 

USED 

BY 

AIRDRG 

AIRFLX 

AIROFF 

AIRP02 

AIR2* 

AIRDAT 

IS 

USED 

BY 

shparf 

ZEROAL 




AIRDAT 

IS 

USED 

BY 

AIRDRG 

AIRFLX 

AIROFF 

AIRP02 

AIR2* 

AIRPOl 

IS 

USED 

BY 

AIRDRG 

AIRPRM 

LSTDY 

LSTUNS 


BDSCOI 

IS 

USED 

BY 

bosoib 





BDSC02 

IS 

USED 

BY 

BDS01S 





BDSC03 

IS 

USED 

BY 

BDS01* 





BESAR 

IS 

USED 

BY 

BE5SAV 

BSINT 

BSJINT 

BSVIWT 


BESDEL 

IS 

USED 

BY 

BESSAV 

BSINT 

8SJJNT 

BSVINT 


BLDGC3 

IS 

USED 

BY 

FETCH 

STORE 




BLDGEO 

IS 

USED 

BY 

BLDCT1 

BLDCTJ 

BLDG 11 

8LDOI2 

BLDGI5 

BLDGEO 

IS 

USED 

BY 

LSTUNS 

MCAFA 




BLDGXO 

IS 

USED 

BY 

FETCH 

RETREV STORE 



BLDOOl 

IS 

USED 

BY 

BLDCT1 

BlDGI 1 

FETCH 

RETREV 

STORE 

BLDG02 

IS 

USED 

BY 

FETCH 

RETREV 

STORE 



BLDGQ3 

IS 

USED 

BY 

BLDCT3 

BLDG 1 3 FETCH 

STORE 


COMO 02 

IS 

USED 

BY 

GGS0I1 

GGS187 




CRPPC1 

IS 

USED 

BY 

CRPPRM 

RETREV 

LSTCTL 

LSTDY 

LSTUNS 

DTETME 

IS 

USED 

BY 

AIRPRM 

BLDOEH 

BLDGI 1 

BLOCI2 

BLDGI3 

DTETME 

IS 

USED 

BY 

WAX0Q1 





DTETHE 

IS 

USED 

BY 

LSTCTl 

LSTDY 

LSTPRM 

LSTUNS 

NEWPOl 

HEADR1 

IS 

USED 

BY 

HEADER 

LSTCTL 

LSTDY 

LSTUNS 

NEWPGi 

HE ADR 2 

IS 

USED 

BY 

HEADER 

LSTCTL 

LSTDY 

LSTUNS 

NEHPGl 

MSQSPL 

IS 

USED 

BY 

LSTSW1 

NOZNHF 

QSPLIN 



INTP01 

IS 

USED 

BY 

INTERS 

INTFRI 

INTPRM 

LSTCTL 

LSTDY 

KCLSAV 

IS 

USED 

BY 

ISOARF 





LSTPOl 

IS 

USED 

tv 

FETCH 

RETREV 

STORE 

LSTDY 

LSTPRM 

LSTP02 

IS 

USED 

BY 

FETCH 

RETREV 

STORE 



LSTRS2 

IS 

USED 

BY 

FETCH 

STORE 




LSTROI 

IS 

USED 

BV 

AIRDRG 

FETCH 

STORE 

INTFRI 

LSTDY 

LSTR02 

IS 

USED 

BY 

FETCH 

STORE 

INTFRI 

LSTDY 

LSTUNS 

LSTSPR 

IS 

USED 

BY 

FETCH 

STORE 




NOZDAT 

IS 

USED 

BY 

NOZCLC 

NO 2 PR M 




NOZSli 

IS 

USED 

BY 

N02FFC 

NO 2 NFC 




OUTPUT 

IS 

USED 

BY 

BDSOl* 





PROVRS 

IS 

USED 

BY 

LSTDY 

LSTUNS 

NEWPOl 



RNCP01 

IS 

USED 

BY 

LSTCTl 

LSTDY 

LSTUNS 

MCAFA 

RUNPRH 

VTXCOH 

IS 

USED 

BY t 

: FETCH 

RETREV 

STORE 

LSTDY 

VTXCAL 

VTXC02 

IS 

USED 

BY : 

FETCH 

RETREV 

STORE 



VlIICl 

IS 

USED 

BY : 

LSTCTl 

LSTDY 

LSTUNS 

VCLGRD 

VEL 002 

WAKPOl 

IS 

USED 

BY : 

FETCH 

RETREV 

STORE 

LSTDY 

LSTUNS 

WAKP02 

IS 

USED 

BY i 

FETCH 

RETREV 

STORE 



WAXR01 

IS 

USED 

BY : 

FETCH 

STORE 

INTFRI 

LSTDY 

LSTUNS 

WAKSP2 

IS 

USED 

BY ■ 

FETCH 

STORE 




WORK 

IS 

USED 

BY : 

F271 

F271NO 

IWTFR3 

LSTCTL 

LSTDY 

WORKl 

IS 

USED 

BY : 

F271 

F271NO 

LSTDY 

LSTUNS 


XINOEX 

IS 

USED 

BY : 

BDSOl* 






CASARF DRAGS* ISOAFt ISOARF LIFT2A 

CASARF DRAG2* ISOAFL ISOARF LIFTS* 

CASARF 0RAG2* ISOAFL ISOARF LIFT2* 

CASARF DRAC2* ISOAFL ISOARF LIFT2* 


BLDOPl FETCH RETREV STORE LSTDY 


LSTWC NEWPG1 NRKOOI WAX001 
BLDOPl HEADER INPTRW INTFR1 INTPRM 

NO I SEX RUNPRH VTXPRH HAXPRM WAKPR1 


LSTUNS LSTWC WAKOOI 
UAK001 


VTXQUT VTXPRH 


VTTBIN 

WAKPRH WAKOOI 


HAKPR1 WAX 081 
LSTUNS NRX801 WAX003 


■ in TSO FOREGROUND HARDCOPY ■■■■ 185 34:35 83112 
DSNAME«TSOOO*I . XREF* . DATA 


LABEL COMMON AREA REFERENCE LISTING 
SUBROUTINE UGES 


5 LABELED COMMON 


5 LSTROi 
: AIROAT 
: AIRDAT 
: AIRP01 
; AIRDAT 
: AIRDAT 
: B DSC 03 
: BDSCOi 
t BESAR 
: BLDOOl 
: BLDC03 
i dtethe 

: DTETHE 
: DTETHE 
J dtethe 
: DTETHE 
i BESAR 
: BESAR 
: BESAR 
s AIRDAT 
; CRPP0I 
5 AIRDAT 
: BLDG03 
; HAKP01 
: WORK 
j WORK 
5 COM802 
: COH082 
j DTETME 
i DTETME 
: INTP01 
: LSTR02 
: WORK 
t DTETME 
: AIRDAT 
j AIRDAT 
j AIRDAT 
: WORK 
: INTP81 
: L5TP81 
: DTETME 
i HSQSPL 
; INTP81 
i LSTP01 
: LSTROI 


AIRP01 AIRDAT AIRCDF 
AIRCDF AIRCDM AIRALZ 
AIRCDF AIRCDM AIRALZ 
DTETHE 

AIRCDF AIRCDM AIRALZ 

AIRCDF AIRCDM AIRALZ 

X INDEX OUTPUT 

BOSC02 

BESDEL 

BLDGEO 

BLDCEO 

BLDOOl BLDGEO 
BLDGEO 

BLDG03 BLDGEO 

BLDCEO 

BESDEL 

BESDEL 

BESDEL 

AIRCDF AIRCDM AIRALZ 

AIRCDF AIRCDM AIRALZ 
BLDOC3 LSTROI LStSPR 
WAXP02 LSTP01 LSTP02 
WORK l 
WORK! 


HEADR1 HEADR2 


AIRCDM AIRALZ 


LSTR02 LSTRS2 WAKRfll WAKSP2 

VTXCOH VTXC02 BLDOOl BLD082 BLDGEO BLDGXO 


RNCP01 CRPP01 INTP01 V1I1CI 

HE ADR 2 DTETHE WAKPQ1 WAXR01 BLDGEO V1I1C1 

LSTROI WORK WORK l VTXCOH 


CRPP01 PROVRS HEADftl 
LSTR02 RMCP01 AIRPOl 
CRPP01 


HE ADR 2 DTETME WAKP01 WAKRfll BLDGEO VII1C1 
LSTROI WORK WORK1 


BLDG01 BLDGEO RNCP01 

PROVRS DTETME HEADR1 HE ADR 2 CRPP01 

DTETME 

NOZDAT 

NOZSli 

NOZSli 

HSQSPL 


Figure 3D: UCAP Labeled Common Area Reference Listing 



Figure 3D: UCAP Labeled Common Area Reference Listing (cont) 
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Figure 6 Blade geometry nomenclature (2-D coords) 
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MFC2140.004 




Figure 7 Blade geometry nomenclature (RXY coord) 
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MFC21 40.005 



I 

Figure 8 Blade geometry nomenclature (RXY coord) 

i 
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MFC21 40.050 
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MFC21 40.007 


▲ Z (Pitch change axis) 



Figure 10 Blade geometry nomenclature (XYZ coords) 


EUnUr"?™":"""!""""" ****«'■*»«*»»** MHMMHNNNMHMMMNNMMMNNNN 

HEADER ( CH"Z i 

UCAP FOR CR-2 INPUT -- R252INPT . DATA 

HX = 0.265, J » 1.^6, TEMP * 56.2, P/PREF = .9627, VTIP = 636 FPS 
6X5 CONFIGURATION, RUN 252.3 - 50/50 POWER SPLIT 
***** FIRST ATTEHPT AT POWER MATCH ***** 
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LSTPARMS(FWD > 

tOOl 1.0 10.0 0.02 6.0 1024. 

C *** PARTI - 0 TO CREATE , 3 TO READ 
L020 1.0 




Figure 11 : • Input E<jho (continued) 



UCAP FOR CR2 INPUT -- R252INPT .DATA 

MX = 0.265, J = 1.46, TEHP = 56.2, P/PREF = .9627, VTIP = 656 FPS 
6X5 CONFIGURATION, RUN 252.3 - 50/50 POWER SPLIT 
***** FIRST ATTEMPT AT POWER HATCH **«*« 


a a 

UJ UJ 

oc cc 
Z> =5 
</></> u 
< < < 
uj uj u 

X E M 


O (J M 
(U IU 

a a a 


S5i 


O N (V O N IT 
. . sO sO LT 

H N 

Iflfr • ■ ■ 


x e e o o o e 


x O O > > O 

x u. r o o z 

XU * (X < < 

*5t\ z 
x m r p o h- y 

X UJ IU I 2 U. < 

X OHKU.C 


l 


I* 
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Figure 12: RU^IPAR^IS output 



a 

M 


O U 
»- I- 
< < 
X X 


t A, u. lu DO 


c h a; o 
i NHL 

: in v 

[ N D. 6 


• m i 
i 5. It w> 

4 C. N 

r m 

UJ M 


IKK 

* S D 

* tf> Ifl J 
X < < < 

UJ UJ u 

= X X * 
O t 
K l i 2 
< LU 

X • • o 

U U H 

Of uj uj 
uiQQO 


H H N O O (NJ M 

IAN O 
N «H 


h II N VI 
- CL UJ 

: x uj uj o 

I UJ -j -i *< 

ihUU J 
• I- Z 2 A 

t < < < 


U '31 
IAO l 
K 'O H I 
ONU.) 


ISel 


X 2 Of 2 

I X 5 < p 
: * Oc uj ec 


H N w ^ ui s 


X O C 3 XX 

U h- 2 K as 

h- M M h" < 

H } u H XX 

X Ifl < X IU 

in #. m U I- Q 
a. iA -J 2 u. X 
< < I 


I V) I 


: A A 


A N O of UJ 
LO.KUI JBH 

as os o *- o x u. 

UUKHhL^ 


a 

ff 


3 

CL 

+* 

3 

O 

t/T 

5 

cc 

< 

a.~ 

a 

cc 

o 


co a 

; ■ 

w 

3 

o> — 

U. ~~ 


1 

1 
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UCAP FOP CP -2 INPUT — R252INPT . DATA 

HX - 0.265, J a 1.46, TEMP - 56.2, P/PREF = .9627, VTIP = 656 FPS 
6X5 CONFIGURATION, RUN 252.5 - 50/50 POWER SPLIT 
NNNM FIRST ATTEMPT AT POWER NATCH ***** 


r~ > 


* JC CL N 

* 5 > j i 

* Zl-U. 

X §iLlL ‘3 

X fl K K K 
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Figure 14: AIRPARKS output 


UCAP FOR CR-2 INPUT — R252INPT . DATA 

MX * 0 . 265, J » 1.46, TEMP *» 56.2, P/PREF * .9627, VTIP « 636 FPS 
6X5 CONFIGURATION, RUN 252.3 - 50/50 POWER SPLIT 
MMMMM FIRST ATTEMPT AT POWER MATCH «■*»* 

FRONT BLADE ANGLE = 98.67 DEG. - MEASURED - CP MATCH = 48.87 


o 

4-» 

o 

T5 

k. 

co 


(O 

a> 

<3 

c 

•o 

w 

O 

cr~ 

u 

O 

■ 


co in 

ui ni 


^ p* to 


o 
o 
v a 

5 

Q. 

5 

o 

O 

LU 

a 

LU 

O 


9 

I 


O O I 

e in i 
no » < 


QQ 


n V) 
iu 
uj a 


oc z 

< o 

LU K 

ac u. 


t • <we ><ANen e o o o po o in • © © ■ o o o 

iiH«owN<r vein c in o o m to in c ph o o cm © in in 
x • <r * • « m m sc in o in © © ** ■ • w © cj cr 

x NcgirviHNiHt-teincgo© o o h o 

* • «8 * © • • • 
x © ©©©©toe©© h o e e 

X CM F I 

X 

X 

X 

« p* p* — . ~ a p* p* — ■ -k 

X pH pH pH pH pH pH pH pH O pH ©pH 

xetototOkTkT^'iTkT in k© n. © ph n to © © m m 

x oooowKitotoe © © © o r*. rv s s <n «< ^ a* 

; - - - -- ™ 
x lu ui V) m 5 5 

x -i -i no mm 

x uu < &. a: os t- t- 

X ui£z t- > fi 2 < < 

X CL < < UJ t-XSL-W 

X K >■ k X U U V) to 

xiuf-OS ►» »jn:m 

xl-zzu oc *-!»»-• 

X u z UJ CL e H IU < o u. 3 

X =) O r IU H 7 UJ » H a. 0 - a. u. Ot- Ot- 

i a<<oiu^:l m d r -i x z s k • cn * to 

X u JHU3U3> V X < UJ W 10 V) M Q W O M 

XD«D</»tOD0CL-X I- A U Q X >■ M < Z -J Z -I 


LO 


3' ’ 

o> 


I 
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AERO PANEL INPUT PA RHS-- FORWARD: LSTPRM LSTPARMS FWD TIHE: 17:00:20 DATE: 1/15/92 UCAP — NASA 


O O O N 

in o o o 
in w o o 


i o a 

: * iu ui 

t * or or 

► * D => 

* 10 (A -J 

> I < < < 

\ U UJ u 

. = r r m 


o o © o 
o o e cm 
in o o o 
to m o o 


e • ■ a 

O CJ H 
) OC UJ UJ 
'UIQQfi 


> o • * o • o e e • • • 

iNn<f'f^ooeeeot 
to CM © H 6 O N 'C rt 

O O O •— 4 LA 
O |H • o o © 

o o o o 


iOHOOWOIfl9><JNOOOON 

© CM (H • CM *40000 
o o o o o cm m o o 


t- O O -J 

Kzza 

< •< < 

or 

mu w < 
tAQOut 
k < < or 


lOOOOOOOOHJrtHJrtHHHHNNNNNNNNNNO'JON 

IOOOOOOOOOOOOOOOOOOCOOOOOOOOOpHHNIM 


W rl N H Z fir 

'ZZZK 3 4 

- H W fl K O k- < 

: or or oc < o </i : 
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, Figure 1 6: LSTPARMS output , 


RUM PARAMETERS FOR MOISE CALC, i MOZPRM NOIZPARM TIME: 17*00:20 DATE: 1/15/92 UCAP — - NASA 


ge 

l if BO 

I I IUUI 
O 1 KK 
A. * P 3 

* i/l ifl 

0 a < < < 

ift umu L) 

v r t r w 

© u j- 

in h i ix 

< IU 

1 r • * o 
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K) o: uj 

* IU o O D 
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:*52P>--<o-ioppxxxo: 

JNNOVIZZ A JU>LXUINNI 


116 


Figure 17: NOIZPARM output 



AXIS. NACELLE AXIAL VEL. *** t VELOOI VELGHADS TIME tMII:MMiSS DATE i fti/OO/YY UAAP 


o tr -i 
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Figure 18 VELGRADS output 
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Figure 18 VELGRADS output (continued) 
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Figure 1 8 VELGRADS output 



VORTEX LOAD PARAMETERS - - FORWARD : VTXPRM VORTPARM FWD TIME; 17:00:20 DATE: 1/15/92 UCAP --- NASA 
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Figure 19: VORTPARM output 



PROPELLER HAKE ANALYSIS INPUT iHAKPRM HAKEPARM TIME .HHiMliSS DATE I W/DD/YY UAAP — NASA 1 
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Figure 20: WAKEPARM output 
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Figure 21: SNTERPRM output 
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Figure 22: AEROEXEC output (GETXMB section 1) 


NUMBER OF INPUT RADII: 22. ROTOR: FORWARD 

STATION RADIUS CHORD TWIST MCA FA 

1 0.100 0.1747 70.004 0.01746 -.00192 




Figure 23: AEROEXEC output (GETXMB section 2) 



Axis of rotation 
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(positive as shown) 
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Figure 24 Blade section definition and definition of advance coordinate, y. 
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Figure 25 AEROEXEC output (getxmb section 3) 
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Figure 25’ AEROEXEC output (getxmb section 3, continued) 
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Figure 26: AEROEXEC output (GETXMB section 4) 
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Figure 27: AEROEXEC output (GETXMB section 5, LESWP output) 
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Figure 28: AEROEXEC output (GETXMB section 7) 
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Figure 29A: AEROEXEC output (CPMAP, front rotor) 
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Figure 29B: AEROEXEC output (CPMAP, rear rotor) 
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Figure 30: AEROEXEC output (HWng surface solver) 


UCAP FOR CR-2 INPUT -- R252INPT . DAT A 

MX * 0. 265, J = 1.46, TEMP * 56.2, P/PREF * .9627, VTIP « 636 FPS 
6X5 CONFIGURATION, RUN 252.3 - 50/50 POWER SPLIT 
mmkxn FIRST ATTEMPT AT POWER MATCH ***** 
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Figure 31 AEROEXEC output (1st camber table) 
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Figure 31 : AEROEXEC output (1st cambor table, continued) 
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Figure 32 : AEROEXEC output (ANSO) 
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Figure 32 : AEROEXEC output (ANSO, continued plus coefficient of sound power) 
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Figure 33: AEROEXEC output (Individual rotor performance summary) 
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Figure 34: AEROEXEC output (vortex calculations: operating conditions) 
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Figure 35: AEROEXEC output (vortex calculations: elemental data) 
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Figure 36: AEROEXEC output (vortex calculations: coefficients from potential calculation) 
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Figure 38: AEROEXEC output (vortex calculations: performance summary) 
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Figure 39: AEROEXEC output (Individual rotor performance summary) 
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Figure 40 : AEROEXEC output (steady interference velocity field calculation) 
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Figure 40 : AEROEXEC output (steady Interference velocity field calculation, continued) 
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Figure 40 AEROEXEC output (steady interference velocity field calculation, continued) 
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Figure 40 : AEROEXEC output (steady Interference velocity field calculation, continued) 
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Figure 41 : AEROEXEC output (PRTPRF: counter-rotation performance summary) 



PERFORMANCE SUHMARY INCLUDING VORTEX LOADING. 
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Figure 41 : AEROEXEC output (PRTPRF: counter-rotation performance summary, continued) 
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Figure 42 : AEROEXEC output (ADJUST) 



POTENTIAL 20-ARY OTHER ROTOR 
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Figure 42 : AEROEXEC output (ADJUST, continued) 
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Figure 42 : AEROEXEC output (ADJUST, continued) 
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Figure 42 : AEROEXEC output (ADJUST, continued) 
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Figure 42 : AEROEXEC output (ADJUST, continued) 
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Figure 42 : AEROEXEC output (ADJUST, continued) 
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Figure 43: AEROEXEC output (CPMAP, rear rotor onto front rotor) 





Figure 44: AEROEXEC output fCPMAP, front rotor onto rear rotor) 
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Figure 45 : AEROEXEC output (steady WMU vector calculation) 


WMU VECTOR! FRONT ROTOR ONTO REAR 
STEADY TERH. 

STATION PHI SICHA1 SIGMA2 SIGHAI/SIGMA2 TWIST 

9 0.5500 99.2903 I .<1133139 1.559253 0.922076 57.1517 
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Figure 45 : AEROEXEC output (steady WMU vector calculation, continued) 
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COMPLEX NUMBERS PRINTED AS (REAL, IMAGINARY) 

fl 0.950 0.029 +29.288 +2.0462 +2.2855 +0.8944 +42.593 

HOAR MU XMDAR DPMI WMU DOWNWASH CMORDWISE VISCOUS HAKE 

1 71 0.506 1.378 ( -10.220, 26.297) t 0.383, -0.102) ( -0.444, -0.010) 0.009 
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Figure 46 : AEROEXEC output (unsteady WMU vector calculation, continued) 
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NOISE CALCULATION OUTPUT : NQ20UT WOIZEXEC EXECCASE TI« 


yCAP FOR CR-2 INPUT — R252IWPT . DATA 

KX • 8.245, J * 1.44* TEMP ■ 54.2, P/PREF - .8427, 
8X5 COFIGlAATIOW , RUN 252.3 - 58/58 POWER SPLIT 
..... FIRST ATTEMPT AT POWER MATCH — 


VTIP 


434 FPS 


FRONT SLADE ANGLE » 48.67 DEG. - MEASURED - CP MATCH ■ 46,67 
REAR SLADE ANGLE ■ <6.54 DEG. - KCASURED - CP MATCH * 46.54 
FRONT l REAR BLADES ASSU8CD IDENTICAL 



NOISE PROGRAM EXECUTION TIME 

CPU TIME AT START « 724.841 SEC 

CPU TIME AT END • 733.448 SEC 

NET TIME - 4.718 SEC 

OBSERVES LOCATIONS 


DIRECTIVITY POINTS AT ALTITUDE OR SIDELINE DISTANCE, 4.120 

DISTANCE FORWARD OF PLANE OF ROTATION (RETARDED) RE FRONT ROTOR 

4.348 4.445 2.573 6.4170 -1.447 -4.111 

DISTANCE FORWARD OF PLANE OF ROTATION (VISUAL) RE FRONT ROTOR 

4.821 2.408 0.6137 - 1.208 -3.114 -4.04< 

RADIATION ANCLE FROM FLIGHT DIRECTION - DEGREES RE FRONT ROTOR 

43.80 52.80 47.20 64.10 183.3 123. ( 

VISUAL RADIATION ANGLE - DEGREES RE FRONT ROTOR 

54.48 44.81 62.43 181.2 117.0 134. 

FRONT ROTOR AZIMUTHAL OBSERVER ANGLE IS 8 8.0 DEGREES 

REAR ROTOR POSITIONS OFFSET FROM FRONT ROTOR OBSERVER POSITIONS BY 


-1.447 

4.112 

-6.773 

FRONT ROTOR 

-3.114 

4.044 

-11.41 

I FRONT ROTOR 

183.3 

123.8 

145.1 

117.8 

134.7 

152.2 


8.2578 ■ FRONT RADIUS 


DISTANCE FORWARD OF PLANE OF ROTATION (RETARDED) RE REAR ROTOR 

7.814 5.274 3.144 8.8588 -0.8444 -3.451 -6.341 

DISTANCE FORWARD OF PLANE OF ROTATION (VISUAL) RE REAR ROTOR 

L547 3.135 1.340 -8.4624 -2.566 -5.548 -11.86 

RADIATION ANGLE FROM FLIGHT DIRECTION - DEGREES RE REAR ROTOR 

41.11 48.23 42.45 81.88 86.78 128.6 143.7 

VISUAL RADIATION ANGLE - DEGREES RE REAR ROTOR 

53.38 42.66 77.45 84.34 112.8 132.2 151.1 


REAR ROTOR AZIMUTHAL OBSERVER ANGLE IS 88.8 DEGREES 


AMBIENT ATMOSPHERIC CONDITIONS 

-•TEMPERATURE 54.15 r , ’‘y...V" 

PRESSURE / SEA LEVEL ■ y 8.8586 - 

* 4* "V 'LOCAL SPEED OF SOUND ■ ~ ~ 1113.34 FEET PER SECOND . > 


Figure 47:NOIZEXEC output (observer position summary) 
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FRONT ROTOR OPERATING CONDITIONS 


8.2^sa FLIGHT MACH NUMBER. 0.5702 TIP ROTATIONAL MACH WV«»CR , 

2.1514 KT / KX, 512.37 BPF , 2.0447 FOOT DIAMETER 


0 .4265 TIP RELATIVE MACH HUHBER , 


rear ROTOR OPERATING CONDITIONS 

- ■ ft, 2451 FLIGHT KACH NUMBER, 1.5733 TIP ROTATIONAL MACH HLMBER, 

2.1*33 MT / MX, 04.33 BPF, 2.0447 FOOT DIAMETER 


0.4314 TIP RELATIVE MACH NUMBER , 


Figure 48:NOIZEXEC output (front and rear rotor operating conditions) 
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MIMT * REA * ROTOR HARKS : WOZOU2 NOIZEXEC EXECCASE TIKE: 17:90:20 DATE: 1/15/12 UCAP HASA < 

"ncAp" " fOR*CR- Z INPUT — R252 INPT .DATA 

!£?- 0^245, J • 1.44, TEMP - 54-2. P/PREF - .1427, VTIP - 454 FPS 
4X5 COMF1 CUR AT 10W, RUN 252.3 - 50/50 POWER SPLIT 

..... first attempt at power match ■ ■■»■ 

FRONT BLADE ANGLE ■ 46.67 DEC. - MEASURED - CP HATCH ■ <6.67 
REAR BLADE ANGLE ■ 46.54 DEG. - MEASURED - CP MATCH ■ <6.54 
FRONT 6 REAR BLADES ASSUMED IDENTICAL 


RESULTS AT OBSERVER POSITION 

FRONT ROTOR 
REAR ROTOR 


VISUAL COORDINATES 
X POSITION ANCLE 
4.02 54.41 

4.55 51.31 


RETARDED COORDINATES 
X POSITION ANCLE 
4.34 43.19 

7.11 41.11 


CALCULATED NOISE FOR HARMONICS IN T>C RANGE OF 1 TIMES BLADE PASSING FREQUENCY 


RADIATION M K 

FREQ COUNTER COUNTER DB 
466.42 2 -1 -45.6 

■ ■ FRONT ROTOR »* 

■ STEADY NOISE ■ 

512.37 I • H.7 

STEADY THICKNESS 12.2 

QUADS UPOLE -200.0 

RADIAL LOADING -200.0 

TOTAL NOISE 101.7 


REAR ROTOR ■■ 

» STEADY NOISE ■ 

414.33 0 1 116.4 

STEADY THICKNESS 
QUADS U POLE 
RADIAL LOADING 
TOTAL STEADY NOISE 

400.71 -1 2 

304.24 -2 3 



FRONT 

ROTOR 


PHASE 

P REAL 

P IMAG 

RMS PA 

0.6 

0.0000 

0. 0000 

0.0000 


101.2 

>0.3741 

1.6143 

1.1332 

142.5 

-0.7775 

0.2445 

0.6150 

9.0 

0.9110 

9.9090 

0.1000 

0.0 

0.0000 

1.0000 

0. 0090 

116.3 

-1.1535 

2.1406 

2.4316 

202.7 

-15.3011 

-4.4173 

14.5117 


CALCULATED NOISE FOR HARMONICS IN THE RANGE OF 2 TI«S BLADE PASSING FREQLCNCY 


RADIATION 

h K 



FRONT 

ROTOR 


FREQ 

COUNTER COUNTER 

D6 

PHASE 

P REAL 

P IMAG 

RMS PA 

1260.71 

3 -1 

-76.3 

301.3 

1.0000 

0.0000 

0.0000 

■ 

■ FRONT ROTOR 






• 

STEADY NOISE ■ 






1164.75 

2 0 

43.4 

346.5 

0.0261 

-O.OOS1 

0.0215 

STEADY THICKNESS 

42.1 

104.6 

-0.0071 

0.1270 

0.0271 


QUADRUPOLE 

-200.0 

0.0 

0.0000 

0.0000 

0.0919 


RADIAL LOADING 

-200.0 

0.0 

0.0000 

8. 0000 

0.0000 


TOTAL NOISE 

43.4 

44.0 

0.0216 

0.0211 

0.1303 

1066.70 

1 1 

146.1 

170.4 

-541.2114 

11.3227 

554.7573 


REAR ROTOR *■ 
• STEADY NOISE ■ 
912.44 0 2 

STEADY THICKNESS 
QUADS UPOLE 
RADIAL LOADING 
TOTAL STEADY NOISE 

614.42 -1 3 

600.57 -2 4 


REAR ROT DR 


DB 

PHASE 

P REAL 

P IMAG 

RKS PA 

102.1 

215.4 

1.2027 

-2.5050 

2.7767 


64.3 

332.2 

0.5447 

-0.1111 

0.4122 

65.6 

323.5 

0.3125 

-0.2312 

0.3667 

-200.0 

9.0 

0.0000 

0 . 1000 

0.0008 

-200.0 

8.0 

0.0000 

o. ooeo 

1 .8100 

12.0 

326.9 

1.4772 

-0.4232 

0.7164 

-134.0 

54.4 

0.0900 

0.0000 

0.8800 

-200.0 

0.0 

0.0000 

0.0000 

a . oooo 



REAR ROTOR 


D 5 

PHASE 

P REAL 

P IMAG 

RKS PA 

47.1 

252.1 

-0.0015 

-0.0046 

0.0050 

151.4 

213.4 

305.2154 

-700.2474 

743.1241 

47.1 

11.7 

-1.0006 

0.1045 

0.0044 

40.1 

137.7 

-0.0143 

0.0141 

0.0221 

-200.0 

0.0 

0.0000 

0.0000 

0.0880 

-200.0 

0.0 

0 . 0000 

0.0000 

0.0000 

42.2 

131.5 

-0.0171 

0.0113 

0.0256 

-123.0 

243.4 

0.0000 

0.0000 

0.0080 

-200.0 

0.0 

0.0000 

0.0000 

0.0100 


Figure 49: NOIZEXEC output (summary of noise results, by directivity) 
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PROGRAM QWMUGEN 


C 

c 

c 

c 

c 

c 

c 

c 


REAL RWMU (20,20) , IWMU(20 ; 20) 

INTEGER TYPE 

#(©) THIS WILL CALCULATE QWMU FOR LSTPARMS 351 TO 550 
TYPE IS 1 FOR SEARS LIFT RESPONSE FOR SINUSOIDAL GUST 
IT IS ASSUMED THE FOLLOWING INFORMATION ARE KNOWN: 

NCP NUMBER OF CHORDWISE PANELS . LSTPARMS 

NSM NUMBER OF RADIAL STATIONS. LSTPARMS 

BD CHORD/DIAMETER RATIO AT NSM STATIONS (BLADEGEO) 

ZAR LOCATION OF THE RADIAL STATIONS (LSTPARMS 101-110) 

J ADVANCE RATIO 

REAL J, BD(10) , ZAR(10) 


WRITE ( * , * ) ' ENTER ADVANCE RATIO ' 

READ ( * , * ) J 

WRITE ( * , * ) ' ENTER NUMBER OF CHORDWISE PANELS' 

READ ( * , * ) NCP 

WRITE (*,*)' ENTER NUMBER OF RADIAL STATIONS ' 

READ ( * , * ) NSM 

WRITE ( * , * ) ' ENTER \NSM, ' RADIAL STATIONS--ASCENDING ORDER' 
READ ( * , * ) (ZAR (I) , 1=1, NSM) 

WRITE ( * , * ) ' ENTER \NSM, ' VALUES OF CHORD / DIAMETER ' 

READ ( * , * ) (BD(I) , 1=1 , NSM) 

WRITE ( * , * ) ' ENTER Q ' 

READ ( * , * ) Q 

WRITE ( * , * ) ' ENTER GUST VELOCITY, DIVIDED BY ROTOR TIP SPEED' 
READ ( * , * ) WZERO 
DO N=1 , NSM 

SIGMA=3 . 14159/J 
DO M=1 , NCP 

XN=-0 . 5+ { FLOAT (M) -0.5) /FLOAT (NCP) 

RWMU (M, N) =WZERO*COS (2 . 0*Q*BD(N) *XN/SIGMA) 

IWMU (M, N) =WZERO*SIN (2 . 0*Q*BD (N) *XN/SIGMA) 

ENDDO 

ENDDO 

WRITE (*, 1010)351, ( (RWMU (M, N) , IWMU(M,N) ,M=1,NCP) , N=1,NSM) 

1010 FORMAT ( 13 , ( T5 , 6 { F9 . 4 , IX ) ) ) 

STOP 

END 


Figure 51 

Listing of sample program to generate QWMU vector for simusoidal gust 
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0<#) CRPBLK. F 1. !.!©(#) 

90/08/15 

16:47:12 


C @ (#) BLOCK DATA FOR COMMON/CRPPOl/ 

C PURPOSE: BLOCK DATA FOR COMMON/CRPPO 1/ 

C MODIFIED TO MAKE DEFAULT VALUE OF CRPTOL 0.01 

BLOCK DATA CRPBLK 


C- -BLOCK DATA FOR CRPPRM 


C- -COMMON BLOCK FOR CRPPARMS 


C 


C 

C 

C 

C 

C 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


REAL AFTHRM, COUNT, 

$ RATIO, ROW, 


CRPBUG , CRPTOL, FWDHRM, 

SPACE, SWITCH 


COMMON/ CRPPO 1 / CRPBUG , 
$ CRPTOL, 
$ RATIO 


SWITCH, 

SPACE, 

COUNT, 

FWDHRM, 

AFTHRM, 

ROW, 


CRPBUG: DEBUG SWITCH. 0: NO DEBUG PRINTOUT 

1 : ECHO INPUT 


SWITCH: CRP/SRP SWITCH. 0: SINGLE ROTATION (DEFAULT) 

1: COUNTER ROTATION 


SPACE: DISTANCE BETWEEN BLADE PITCH CHANGE AXES, DIVIDED BY 

FRONT ROTOR DIAMETER. 

RATIO: REAR DIAMETER DIVIDED BY FORWARD DIAMETER 

COUNT: MAXIMUM NUMBER OF FRONT/REAR ITERATIONS. DEFAULT IS 0, 

WHICH IS ONE TRIP THROUGH. 

CRPTOL: MINIMUM VALUE OF THE CONVERGENCE MONITOR REQUIRED TO 

START ANOTHER ITERATION. IF THIS IS <= 0, IT WILL BE 
SET TO 1 . OE+4 

ROW: ROW SELECTOR SWITCH: SET INTERNALLY ONLY i 

<1.5: FRONT ROW 

>1.5: REAR ROW 


FWDHRM: HIGHEST FORWARD ROW WAKE HARMONIC TO USE FOR EXCITATION 

OF REAR ROTOR. NOTE THAT THE FORWARD ROTOR WAKE IS 
STEADY WHEN VIEWED FROM A POINT MOVING WITH THE FRONT 
ROTOR. THE HARMONICS ARE HARMONICS IN SPACE ONLY. 

THOSE HARMONICS > 0 WILL APPEAR UNSTEADY TO THE REAR 
ROTOR IN THE SAME WAY BUMPS ON A ROAD CAUSE TIME VARYING 
FORCES ON AN AUTOMOBILE MOVING DOWN THE ROAD. 

AFTHRM: HIGHEST REAR ROW POTENTIAL FIELD HARMONIC TO USE FOR 

EXCIT AT I ON OF THE FORWARD ROTOR. NOTE THAT THE REAR 
ROTOR POTENTIAL FIELD IS, LIKE THE WAKE, STEADY ONLY 
WHEN VIEWED FROM A POINT MOVING WITH THE REAR ROTOR. IT 
IS ALSO NOT HOMOGENOUS. 

DATA AFTHRM/0 .0/, COUNT/0.0/, CRPBUG/1.0/, 

$ CRPTOL/1 . OE-2/, FWDHRM/ 0.0/, SPACE/-1.0/, 

$ SWITCH/ 0.0/ 

END 


Figure 52 

LiBting of CRPBLK with changed default valve for CRP tolerance CRPTOL 
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